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1.  Introduction 

Absorption  by  the  CO2  and  H2O  molecules  dominate  the  infrared  transmission 
spectrum  of  the  atmosphere.  Knowledge  resulting  from  the  research  carried  out  under 
this  contract  has  contributed  to  our  understanding  of  several  different  atmospheric  ^ 
phenomenon.  The  global  heat  balance  of  the  atmosphere  is  dominated  by  absorption 
and  emission  of  H2O  and  CO2,  the  so  call  "greenhouse"  effect.  A  detailed  knowledge  of 
the  spectra  of  the  major  atmosphere  molecules,  also  makes  it  possible  to  detect  the 
spectral  contribution  of  minor  atmospheric  species.  An  additional  application  of  this 
knowledge  is  detecting  high  temperature  targets  through  the  atmosphere.  As  part  of  this 
contractual  effort,  over  12,000  molecular  transitions  frequencies  of  CO2  and  H2O  were 
measured  and  submitted  for  inclusion  in  the  PL/GPOS  HITRAN  molecular  database. 
The  HITRAN  database  is  used  by  the  U.S.  Air  Force’s  atmospheric  transmission 
computer  codes  and  by  the  atmospheric  community. 


2.  Instrumentation  Description 

The  experimental  measurements  were  made  using  the  Air  Force  Geophysics 
Laboratory  high  resolution  Michelson  interferometer  (1,  2)  and  high  temperature 
absorption  cell  (3).  The  high  resolution  interferometer  is  a  step-and-hold  Michelson 
interferometer  employing  cat’s  eye  retro-reflectors  for  both  the  moving  and  fixed  mirrors. 
A  maximum  optical  path  of  114  cm  was  used,  resulting  in  an  apodized  resolution  of 
0.009  cm The  1.75  m  path  stainless-steel  high  temperature  absorption  cell  is  housed  in 
a  commercial  electric  furnace  which  can  be  heated  to  a  maximum  of  1000  K.  Instrument 
control,  data  acquisition,  and  data  processing  were  performed  using  personal  computers 
running  custom  software.  The  spectrometer  and  absorption  cell  has  been  described 
previously  (1,  2).  During  the  course  of  this  contract,  the  computer  controlled  servo 
system  that  controls  the  moving  cat’s  eye  retro-reflector  was  updated  with  additional 
hardware  and  software. 

An  additional  real  time  data  acquisition  board  was  added  to  the  computer  that 
controls  the  high  resolution  spectrometer.  This  board  was  actually  a  complete  real  time 
computer  with  its  own  microprocessor,  analog  to  digital,  and  digital  to  analog  converters. 
It  was  necessary  to  modify  the  control  electronics  to  interface  this  additional  coniputer  to 
the  rest  of  the  servo  ^stem.  Software  was  written  to  keep  both  processors  working 
together  to  control  the  motion  of  the  interferometer  moveable  cat’s  eye  and  record  the 
interferogram  signals.  Software  compatibility  is  much  more  difficult  to  achieve  for  real 
time  control  applications  that  have  to  interact  closely  with  the  hardware,  than  it  is  for  the 
more  general  computational  computer  applications.  These  enhancements  increased  the 
precision,  dependability,  and  diagnostic  capability  of  the  servo  system  that  controls  the 
high  resolution  spectrometer. 
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A  review  of  how  the  computer  controlled  servo  that  controls  the  motion  of  the 
movable  cat’s  eye  is  helpful  in  understanding  the  enhancements  that  were  made  to  the 
control  ^stem.  The  interferometer  is  a  step-and-hold  instrument  that  uses  cat’s  eye 
retro-reflectors.  The  maximum  optical  path  length  the  interferometer  can  achieve  is 
more  than  1  meter.  To  achieve  these  large  changes  in  optical  path  length  while 
maintaining  the  necessary  optical  path  length  positional  accuracy  needed  for 
interferometric  measurements,  the  interferometer  uses  three  sources  of  motion.  These 
three  sources  are:  1)  a  mechanical  lead  screw,  similar  to  that  used  in  a  metal  working 
lath,  is  used  to  provide  large  motions,  2)  a  linear  electric  motor  is  used  to  provide 
motions  of  the  order  of  1  centimeter,  and  3)  a  piezoelectric  stack  mounted  on  the 
secondary  of  one  of  the  cat’s  eyes  is  used  to  provide  very  fine  optical  control.  A  LVDT 
(linear  variable  differential  transformer)  is  used  to  sense  the  position  of  the  linear  motor 
and  cause  the  lead  screw  to  turn,  maintaining  the  linear  motor  in  the  center  of  its 
working  range.  Fine  positional  information  is  obtained  from  detecting  the  interference 
fringes  from  a  HeNe  laser. 

Before  the  update  to  the  control  system,  the  lead  screw  could  be  moving  while 
data  was  being  acquired  and  the  linear  motor  was  stepping.  The  linear  motor  and  piezo 
servos  were  not  able  to  completely  compensate  for  the  vibrations  caused  by  the  rotating 
screw.  With  the  updated  software,  the  screw  is  not  allowed  to  move  independently  of 
the  linear  motor  and  piezo  stack,  but  instead  is  maintained  in  one  position  until  a  block 
of  interferogram  data  has  been  acquired  and  written  to  the  hard  disk.  Then  the  screw  is 
activated  and  the  linear  motor  is  brought  back  to  the  center  of  its  working  range.  With 
the  new  servo  configuration,  lead  screw  motion  does  not  perturb  the  data  taking  or  the 
stepping  of  the  interferometer.  Problems  of  backlash  and  gear  slap  in  this  lead  screw 
and  its  drive  train  have  also  been  minimized  by  using  the  computer  to  optimize  the 
voltage  that  drives  the  lead  screw  motor.  In  particular,  it  has  proven  veiy  helpful  to  give 
the  lead  screw  motor  a  small  reverse  pulse  at  the  end  of  lead  screw  motion  to  maintain 
the  gear  train  under  load. 

Other  additions  to  the  control  system  of  the  interferometer  were  enhanced 
diagnostics  and  self-correction  for  some  types  of  mirror  motion  problems.  The  new 
system,  using  two  data  acquisition  computers,  makes  possible  continuous  monitoring  of 
the  interferometer  optical  path  difference.  With  the  old  tystem,  the  optical  path 
difference  could  not  be  monitored  at  the  same  time  that  the  infrared  data  was  being 
acquired.  The  optical  path  difference  needs  to  be  continually  monitored  since  external 
vibration  can  occur  at  any  time.  When  the  computer  senses  improper  holding  of  the 
optical  path  difference,  diagnostic  information  is  recorded  to  the  hard  disk  and  the 
affected  data  point  is  remeasured  after  a  short  delay  to  allow  any  external  vibrations  to 
subside. 

The  data  collection  and  reduction  software  is  organized  as  a  number  of  small 
programs  instead  of  one  large  program.  These  programs  communicate  with  each  other 
through  the  headers  at  the  start  of  each  data  file.  The  structure  of  these  headers  was 
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changed  from  a  fixed  number  of  parameters  to  a  variable  number  of  parameters  to 
facilitate  future  enhancements  to  the  software. 


3.  Molecular  Transition  Identification  Techniques 

Probably  the  biggest  problem  in  working  with  the  high  density  spectra,  that  result 
from  high  temperature  gases,  is  making  correct  line  identifications.  In  making  these 
identifications,  a  combination  of  a  variety  of  different  techniques  proved  more  useful 
than  any  single  technique.  A  computer  generated  Loomis-Wood  diagram  proved  the 
most  helpful  in  providing  initial  identifications.  Verification  of  identifications  was 
accomplished  most  effectively  with  combination-difference  methods.  For  the  CO2 
molecule,  some  success  was  obtained  using  an  automated  search  through  the  rotation- 
vibration  parameter  space.  Computer  software  was  developed  to  implement  each  of 
these  identification  techniques. 

3.1.  Looms-Wood  Diagrams 

A  Looms-Wood  diagram  is  used  to  display  spectral  information  graphically  in  a 
manner  such  that  it  is  easy  to  pick  out  the  lines  belonging  to  a  particular  rotation- 
vibration  band  or  sub-band  in  the  presence  of  many  overlapping  bands.  A  version  of  the 
software  needed  to  display  Looms-Wood  diagrams  had  been  developed  previously  (4), 
but  during  this  contract  was  adapted  to  take  advantage  of  the  many  more  colors  that  are 
possible  to  display  using  modem  personal  computers. 

The  first  step  in  the  application  of  the  Loomis-Wood  technique  is  to  compile  a  set 
of  first-guess,  self-consistent  spectral  line  positions  for  the  branch  or  sub-band  whose 
lines  are  to  be  identified.  For  the  original  Looms-Wood  diagram  (5),  the  first-guess  line 
positions  had  constant  spacing.  A  computer  generated  map  is  then  constructed  in  the 
following  way:  1)  equal  wavenumber  segments  of  the  transmittance  spectra  being  studied 
are  selected,  2)  each  segment  is  centered  on  the  predicted  position  of  one  of  the  spectral 
lines  in  the  sub  band  being  analyzed,  3)  the  segments  are  then  stacked  one  above  the 
other.  If  the  predicted  and  observed  spectral  line  positions  agree,  then  the  observed 
spectral  line  will  be  centered  in  the  figure  forming  a  straight  line.  The  set  of  spectr^ 
lines  forming  the  sub-band  will  be  on  top  of  one  another  and  form  an  easily  discernible 
straight  line  pattern. 

If  the  first-guess  transitions  are  in  error,  an  irregular  or  curved  pattern  will  be 
observable.  The  errors  from  most  calculations  follow  some  sort  of  pattern.  For  example, 
if  rotational  constants  are  used  to  determine  the  first-guess  line  positions,  and  the 
quadratic  dependence  on  the  rotational  quantum  number  J  is  in  error,  then  a  quadratic 
deviation,  as  a  function  of  J,  from  the  straight  line  pattern  will  be  observed.  The 
Loomis- Wood  technique  makes  it  easy  to  discover  branch  or  sub-band  spectral 
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lines,  even  in  dense  spectra,  and  to  find  and  correct  differences  between  observed  and 
predicted  line  positions,  if  a  reasonable  set  of  initial  line  positions  are  available. 

The  original  Looms- Wood  diagram  (5)  was  only  useful  for  linear  molecules,  but  it 
has  been  possible  to  extend  this  technique  to  asymmetrical  molecules  like  H2O.  This  is 
done  by  finding  sequences  of  lines  where  the  intensities  and  positions  follow  regular 
patterns.  These  sequences  are  not  nearly  as  long  as  they  are  for  a  linear  molecule  like 
CO2,  but  especially  in  the  case  of  high  temperature  spectra,  they  are  clearly  visible.  The 
sof^are  that  was  developed  to  make  line  assignments  and  display  Looms-Wood 
diagrams  is  interactive.  A  crosshair  is  moved  across  the  diagram  and  used  to  select 
spectral  lines.  The  software  then  automatically  determines  the  line  centers  of  the 
e^qperimental  feature  and  creates  a  new  HITRAN  format  data  file  with  the  updated  line 
parameters. 

32.  Verifying  Line  Assignments 

A  procedure  was  developed  to  check  for  self-consistency  among  the  measurements 
made  on  different  rotation-vibration  bands,  including  bands  from  different  spectral 
regions.  One  can  think  of  following  paths  through  the  energy  level  space,  where  the 
allowed  steps  in  this  path  are  the  observed  rotation-vibrational  transitions.  Starting  at  a 
given  energy  level,  the  energy  of  the  next  level  in  the  path  is  obtained  by  simply  adding 
or  subtracting  (depending  on  the  direction  of  the  path)  the  energy  of  the  transition.  For 
a  closed  loop,  the  ending  level  is  the  same  as  the  beginning  level  and  the  net  energy 
difference  after  traversing  the  path  should  be  zero.  If  experimentally  measured 
transitions  are  used,  the  net  energy  of  the  path  will  not  be  quite  zero,  but  instead  will  be 
a  measure  of  the  relative  accuracy  of  the  measurements.  Calculating  the  net  energy 
around  a  closed  loop  of  energy  levels  is  also  a  veiy  sensitive  technique  for  detecting  line 
assignment  errors.  If  an  error  in  the  line  assignment  has  been  made,  a  path  that  is 
assumed  to  form  a  closed  loop  will  actually  not  be  closed,  and  the  resulting  net  energy  of 
the  path  will  be  much  greater  than  that  expected  from  random  experimental  fluctuations. 
Since  finding  these  closed  loops  manually  proved  very  tedious  and  time  consuming, 
software  was  written  to  automaticjally  generate  these  closed  loops  fi'om  data  files 
containing  the  observed  molecular  transitions. 

Finding  appropriate  sums  and  differences  of  observed  transitions  was  also  used  to 
predicted  the  position  of  molecular  transitions  that  have  not  yet  been  identified.  The 
position  of  a  spectral  transition  can  be  predicted  if  a  path  of  connecting  observed 
transitions  can  be  found  from  the  lower  state  to  the  upper  state.  Software  was  written  to 
automatically  make  an  exhaustive  search  for  all  such  paths  less  than  a  maximum  number 
of  steps  (typically  6  to  7  transitions).  Allowing  longer  paths  would  have  increased  the 
number  of  paths  that  could  have  been  identified,  but  (due  to  the  accumulation  of  errors) 
would  have  resulted  in  less  accurate  predictions  of  molecular  transitions.  Paths  of  length 
6  to  7  steps  proved  to  be  a  good  compromise.  When  several  alternate  paths  were  found, 
a  weighted  average  of  the  transition  frequency  predicated  by  the  different  paths  was 
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calculated.  The  weight  for  each  path  was  taken  as  the  inverse  sum  of  squares  of  the 
individual  uncertainties  of  the  transitions  making  up  the  path. 

33.  Automated  Identification 

Another  identification  technique,  which  was  useful  for  finding  the  rotation- 
vibration  bands  of  CO2,  was  an  automated  search  of  parameter  space.  It  used  an 
iterative  least-squares-fits  method  to  search  for  the  set  of  molecular  constants  that  best 
fit  the  experimental  data.  A  difficulty  with  this  approach  is  that  the  error  function  being 
minimized  has  a  large  number  of  local  minimums,  making  it  difficult  to  find  the  absolute 
minimum.  A  successful  approach  was  to  start  the  search  for  a  minimum  from  many 
different  positions,  finding  many  local  minima  but  only  recording  the  molecular 
parameters  for  the  fits  that  appeared  the  most  promising.  In  this  way,  a  comprehensive 
search  of  the  whole  parameter  space  was  carried  out.  The  most  promising  fits  were  then 
manually  examined  in  detail  using  the  usual  interactive  techniques.  Perfo^ng  large 
number  of  trial  fits  consumes  a  great  deal  of  computer  time.  However,  with  the  cost  of 
computing  hardware  coming  down,  and  with  the  rapid  increase  in  the  computational 
power  of  the  new  desktop  computers,  this  technique  appears  to  have  a  promising  future. 


4.  Molecular  Modeling 

In  addition  to  the  experimental  measurements  of  COj  and  HjO,  extensive 
molecular  modeling  was  also  performed  as  part  of  this  contractual  effort.  Although 
there  are  a  number  of  possible  modeling  techniques,  (DND)  direct  numerical 
diagonalization  (6)  has  proved  the  most  successful  in  calculating  the  extensive  sets  of 
energy  levels  need  to  support  the  identification  of  molecular  transitions  of  high 
temperature  gases.  Calculations  using  the  DND  technique  can  also  be  used  to  ^culate 
line  parameters  for  molecular  transitions  that  are  very  difficult  to  observe  experimentally. 


Work  was  performed  on  CO2  pentad  band  parameters  in  conjunction  with 
researchers  from  other  institutions.  This  work,  for  the  first  time,  validates  DND  band 
parameter  calculations  for  components  of  polyads  where  not  all  of  the  bands  have  been 
observed.  A  new  CO2  HTTEMP  database  was  generated  which  incorporates  these  new 
measurements  and  level  identifications  at  the  end  of  each  line  record.  The  extra  level 
identifications  allows  for  the  vibrational  identification  of  lines  whose  levels  are  not 
defined  in  HITRAN. 

The  computer  modeling  software  that  is  used  to  fit  molecular  potential  surfaces  to 
the  experimental  data  was  upgraded  by  symmetry-adapting  the  fitting  procedure  to  utilize 
the  C2V  symmetry  of  triatomic  molecules  such  as  H2O,  CO2,  and  O3.  This  makes 
rovibrational  level  identification  easier  and  more  certain  and  reduces  computer  resource 
requirements.  The  updated  software  was  then  used  to  find  a  stable  quartic  potential 
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surface  for  water.  A  new  hot  water  database  was  also  generated.  This  database  achieves 
better  convergence  for  high  combined  vibrational  and  rotational  states. 

An  HjS  spectroscopic  database  has  been  generated  using  molecular  surfaces  obtained 
from  recent  work  of  Kozin  and  Jensen  (7)  (potential  surface)  and  re-adjustment  of  the 
dipole-moment  surfece  published  by  Carter  et  al.  (8).  Comparisons  with  the  new 
HTTRAN  database  for  H2S  show  that  there  are  line  identification  problems  and  possibly 
line  intensity  problems  with  the  HITRAN  92  databzise. 


5.  CO2  Measurements 

The  details  of  the  COj  measurement  program  are  described  in  journal  articles 
included  in  the  Appendix  and  will  not  be  duplicated  here. 


6.  H2O  Measurements 

A  20  torr  sample  of  natural  isotopic  abundance  H2O  was  heated  to  lOOOK  in  the 
high  temperature  absorption  cell  and  observed  with  the  high  resolution  interferometer 
spectrometer.  A  background  spectrum  was  first  collected  with  the  sample  cell  evacuated. 
The  background  spectrum  was  collected  at  the  same  resolution  as  the  sample  spectrum. 
For  phase  correction,  a  short  two-sided  spectrum  of  123  cm'^  resolution  was  collected. 
The  sample  spectrum  was  collected  over  a  12-hour  period,  of  which  approximately  8 
hours  was  actual  data  collection  time.  The  remaining  time  was  mirror  step  and  settle- 
down  time.  The  resulting  signal-to-noise  ratio  (SNR)  was  better  than  100  to  1  in  the 
center  of  the  spectral  region  studied,  near  1050  wavenumbers  (see  Figure  1).  At  the 
high  wavenumber  end  of  the  spectrum  the  SNR  is  lower,  as  a  result  of  the  optical  filter 
cutoff.  At  the  low  wavenumber  end  of  the  spectrum  the  system  sensitivity  decreases  as  a 
result  of  reduced  transmission  by  the  numerous  zinc  selenide  windows  in  the  optical 
path. 


6.1.  Data  Processing 

A  line-finding  computer  program  was  used  to  determine  the  spectral  line  positions 
in  the  transmittance  spectrum.  The  program  takes  into  account  the  spectral  line 
^mmetiy  and  the  effect  of  blending  by  other  spectral  lines.  A  total  of  approximately 
1600  spectral  line  positions  were  determined  from  the  data  shown  in  Figure  1.  The  line 
position  uncertainties  were  determined  by  the  magnitude  of  the  peak  absorptance  of  the 
spectral  line,  the  width  of  the  line,  the  signal-to-noise  ratio  in  the  region  around  the 
spectral  line,  and  the  degree  to  which  the  line  was  blended  with  other  lines. 

The  sample  water  vapor  concentration  in  the  absorption  cell  was  chosen  to 
achieve  large  absorptance  in  the  weaker  spectral  lines.  This  resulted  in  the  stronger  lines 
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Compressed  plot  of  the  measured  spectrum  of  a  20  Torr  HjO  sample  heated  to  1000K. 


being  saturated  near  line  center.  Since  these  strong  lines  have  been  studied  previously 
(9)  and  are  well  known,  the  measured  positions  of  these  lines  are  not  reported  here. 
However,  these  strong  lines  were  used  to  calibrate  the  wavelength  scale  of  the 
spectrometer,  because  these  lines  were  present  in  the  background  spectrum  of  the 
absorption  cell.  Even  after  being  evacuated,  the  cell  contained  a  small  residual  amount 
of  water  vapor,  and  the  strongest  lines  were  present  in  the  background  cell  spectra. 
Forty-nine  background  spectral  lines,  spread  across  the  observed  spectral  region  from 
1225  to  1398  cm'^  were  used  for  calibration.  The  standard  deviation  of  the  calibration 
residuals  was  3.9x10'^  cm'^ 

The  1600  measured  spectral  line  positions,  were  combined  with  observed 
transitions  from  other  studies,  to  form  a  self-consistent  set  of  transitions  which  could 
then  be  used  to  determine  energy  levels  and  transition  identifications.  The  additional 
transitions  included  162  lines  calculated  from  energy  levels  reported  by  Toth  (10),  and  70 
microwave  spectral  lines  measured  by  Pearson  et  al.  (11).  The  most  difficult  step  in  the 
data  processing  was  correctly  identi^ng  the  molecular  transitions.  Of  the  1600 
measured  spectral  line  positions,  a  total  of  921  ro-vibrational  transition  identifications 
were  made  with  certainty.  Many  additional  transition  assignments  were  made,  but  these 
could  not  be  confirmed  by  the  methods  outlined  below. 

A  number  of  different  techniques  were  used  to  make  the  transition  assignments. 
For  very  well  know  H2O  spectral  lines  it  was  possible  to  unambiguously  pick  the 
observed  transition  closest  to  the  position  given  in  the  HITRAN  data  base  (9).  But  for 
most  of  the  transitions  discernible  at  high  temperatures,  more  sophisticated  techniques 
were  required  to  unambiguously  make  the  identification.  An  enhanced  Loomis- Wood 
diagram-like  technique,  combined  with  DND  results  (6)  proved  extremely  valuable  in 
making  line  assignments. 

An  automated  combination  difference  approach  also  proved  very  useful  in 
verifying  correct  identifications.  The  technique  was  also  used  to  identify  additional 
transitions,  once  a  number  of  the  transitions  had  been  successfully  identified.  The 
identification  was  made  in  an  iterative  manner.  First,  a  number  of  lines  were  identified. 
Knowledge  of  the  energy  levels  derived  from  these  identified  lines  was  used  to  identify 
additional  lines.  The  process  could  then  be  repeated,  allowing  the  further  identification 
of  additional  lines. 

In  past  studies  of  carbon  dioxide  (4,  12),  an  enhanced  version  of  the  original 
Loomis-Wood  diagram  technique  (5)  was  used  to  identify  spectral  lines  in  P,  Q,  and  R 
branches.  At  high  temperatures  the  spectrum  is  crowded  with  many  more  spectral  lines 
than  at  room  temperature.  The  Loomis-Wood  technique  was  veiy  useful  in  identifying 
spectral  lines  of  a  particular  band.  It  has  been  possible  to  extend  our  modified  Loomis- 
Wood  diagram  technique  to  the  asymmetrical  molecule  HjO.  The  rotational  energy 
states  of  HjO  are  generally  enumerated  using  the  quantum  number  J  for  the  total 
angular  momentum,  and  K,  and  for  components  of  the  angular  momentum.  In  order 
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for  the  Looms-Wood  diagram  technique  to  be  effective,  sequences  of  transitions  need  to 
be  calculated  where  the  line  intensities  and  positions  follow  regular  patterns.  The  HjO 
sub-bands  for  a  given  and  constant  AK^  and  AJ  are  examples  of  these  sequences. 

These  sequences  are  not  as  long  as  for  a  linear  molecule  like  CO2,  but  especially  in  the 
case  of  high  temperature  spectra,  they  are  clearly  discernible  using  the  Loomis-Wood 
technique. 

The  first-guess,  self-consistent  spectral  line  positions  needed  to  create  a  Loomis- 
Wood  diagrzims  were  obtained  either  from  the  HITRAN  data  base  (9),  or  from  an 
application  of  the  direct  numerical  diagonalization  (DND)  calculation  technique.  The 
DND  calculation  for  H2O  used  essentially  the  same  method  (6)  as  was  previously  used 
for  CO2.  A  47-parameter  potential  surface,  with  up  to  hextic  powers  in  internal 
coordinates,  was  generated  with  a  nonlinear  least-squares  fit  to  observed  data.  The 
residuals  of  the  fit  result  in  an  RMS  of  0.03  cm'^  for  the  263  rovibrational  energy  levels 
used  as  input  to  the  calculation. 

62.  Results 

From  the  observed  H2O  molecular  transitions,  a  set  of  self-consistent  transitions 
was  calculated  using  an  iterative  approach.  In  addition  to  using  the  observed  transitions 
from  this  work,  selected  transitions  from  a  room  temperature  H2O  measurement  by  Toth 
(10)  and  microwave  transitions  from  Pearson  (11)  were  also  included.  Combining 
information  from  different  regions  of  the  electromagnetic  spectrum  improves  the 
accuracy  of  the  calculation.  This  happens  not  only  because  of  the  increase  in  accuracy 
that  additional  measurements  provide,  but  also  by  reducing  the  correlation  in  the 
calculation.  For  each  iteration,  the  position  of  a  transition  was  calculated  from  a 
weighted  average  of  the  positions  predicted  using  combination  differences  of  the  other 
transitions.  In  the  next  iteration,  the  transitions  were  replaced  with  the  positions 
calculated  in  the  preceding  iteration.  This  process  was  iterated  until  convergence  was 
achieved. 

A  total  of  921  lines  were  identified  in  the  experimental  spectrum.  Actually, 
considerably  more  than  the  921  transitions  reported  were  observed,  but  it  was  only  these 
lines  that  could  be  verified  using  the  combination  difference  technique  explained  earlier. 
Positional  uncertainties  were  estimated  from  the  quantity  of  random  noise  in  the 
experimental  spectral,  line  strengths,  and  blending  of  unresolved  lines.  The  amount  of 
line  blending  from  unresolved  lines  was  estimated  from  asymmetry  in  the  line  profiles 
and  inconsistent  line  width  for  a  given  intensity. 

When  the  H2O  line  positions  determined  in  this  work  were  compared  with  the 
line  positions  on  the  HITRAN  92  database  (9),  discrepancies  of  up  to  0.8  cm'^  were 
found  (Figure  2)  for  some  transitions.  These  large  discrepancies  probably  arose  from  the 
inclusion  of  data  from  misidentified  H2O  lines  on  the  HITRAN  92  database.  In  the 
bottom  panel  of  Figure  2,  the  vertical  scale  has  been  expanded  by  a  factor  of  10.  It  can 
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Figure  2.  Coruparison  of  observed  H2O  line  positions  with  line  positions  from  the 
HTTRAN  92  database.  The  vertical  scale  has  been  expanded  by  a  factor  of  10  in  the 
bottom  panel. 


be  seen  that  there  are  a  number  of  additional  lines  where  the  difference  is  smaller  than 
0.1  cm'^,  but  still  large  compared  to  the  experimental  uncertainty  of  about  0.0004  cm  . 
The  HITRAN  database  is  an  evolving  source  of  data  on  atmospheric  molecules.  Hence, 
H2O  lines  from  the  present  work  have  been  incorporated  into  the  HITRAN  96  database. 
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Abstract — More  than  6500  P,  Q  and  R  branch  spectral  lines  in  the  region  580-940  cm"‘  have 
been  identified  and  assigned  to  43  vibration-rotation  bands  of  the  principal  isotope  of  carbon 
dioxide,  The  position  of  each  line  was  measured  and  an  uncertainty  in  the  position 

of  each  line  was  estimated.  The  interferometric  data  used  in  this  study  were  obtained  from 
naturally  occurring  samples  of  carbon  dioxide  heated  to  800  K,  to  yield  data  on  high 
quantum-number  rotational,  vibrational  and  vibrational  angular  momentum  states.  Measure¬ 
ments  were  made  for  three  different  sample  pressures:  5,  15  and  40torr.  Measurements  at  15 
and  40  torr  permitted  the  observation  of  weak,  high  quantum-number  vibrational  transitions 
normally  not  observed  at  room  temperature  and  high-quantum-number  rotational  transitions 
within  the  vibrational  transitions.  A  maximum  unapodized  instrumental  resolution  of 
0.0044  cm“^  was  used.  A  weighted  linear  least-squares  procedure  was  applied  to  the  spectral 
lines  of  each  vibrational  transition  to  estimate  rotational-vibrational  molecular  constants 
B.,,  and  which  are  reported  here. 


INTRODUCTION 

The  objective  of  the  current  study  is  the  determination  of  spectral  line  positions  of  carbon  dioxide. 
The  goal  is  to  collect  spectral  data  with  the  largest  possible  number  of  observed  spectral  lines 
corresponding  to  high  quantum-number  rotational,  vibrational  and  vibrational  angular  momentum 
transitions.  The  set  of  assigned  and  measured  line  positions  resulting  from  the  data  analyses  may 
then  be  combined  with  data  from  other  sources  and  used  with  available  theoretical  and 
computational  methods'  in  the  determination  of  a  much  larger  set  of  line  positions  for  inclusion 
in  the  HITRAN^-^  and  HITEMP^^  data  bases.  To  study  spectral  lines  corresponding  to  high 
quantum  numbers,  the  gas  sample  was  heated  to  high  temperature.  This  populates  high 
quantum-number  vibrational  and  rotational  energy  levels  and  yields  many  more  bands  and  more 
spectral  lines  in  each  band  than  has  been  available  before.^'^ 

EXPERIMENTAL  DETAILS 

The  instrumental  arrangement  has  been  discussed  in  detail  elsewhere'^  and  will  be  described  only 
briefly  here.  The  interferometric  measurements  were  made  with  a  step-and-hold  Michelson 
interferometer  employing  cat’s-eye  moving  and  fixed  mirrors.  The  maximum  stepping-mirror 
displacement  for  this  study  was  57  cm  corresponding  to  an  unapodized  resolution  of  0.0044  cm“'. 
The  output  beams  of  the  cat’s-eye  mirrors  are  laterally  displaced  making  both  accessible  and 
allowing  the  use  of  two  detectors.  Because  the  two  beams  are  complementary,  the  detector  signals 
are  differenced  to  yield  the  interferometric  signal.  This  arrangement  reduces  the  effect  of  unwanted 
external  disturbances  such  as  60  Hz  interference,  source  fluctuations  and  instrumental  drifts.  The 
arrangement  provides  a  very  stable  system  over  the  8-15  h  period  necessary  for  the  collection  of 
a  single  high-resolution  interferogram.  Both  liquid-helium  cooled  copper-doped  germanium 
detectors,  used  for  operation  in  the  15  jum  spectral  region,  are  housed  in  the  same  dewar  to  keep 
their  experimental  conditions  the  same.  An  interference  filter  was  used  to  limit  the  spectral  range 
from  580  to  940  cm 'L  The  interferometer,  HeNe  reference  laser,  signal  detectop,  optical  filter 
wheel  as  well  as  all  of  the  infrared  and  laser  transfer  optics  between  the  absorption  cell  and  the 
interferometer  are  housed  in  the  same  chamber  which  is  evacuated  during  operation. 
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A  1  75  m  stainless-steel  gas-sample  cell  was  used  in  this  study.  It  was  housed  in  a  commercial 
electric  furnace  which  can  be  heated  to  at  least  800  K.  This  cell  (used  in  a  triple  pass  Pfund 
configuration)  and  furnace  have  been  described  previously.''*  For  the  current  set  of  measurements 
a  single  pass  through  the  cell  was  used.  This  eliminated  optical  misalignment  problems  which  result 
from°changes  in  the  cell  length  with  changing  temperature,  approx.  1  cm  from  room  temperature 
to  800  K.  Only  the  central  windows  are  at  room  temperature.  However,  since  only  measurements 
of  spectral  line  positions,  and  not  line  strengths  or  widths,  were  the  object  of  this  study,  no  attempt 
was  made  to  characterize  these  temperature  gradients.  This  instrumental  difficulty  has  been 
considered  in  a  separate  study.'^ 

The  source  chamber  contains  a  Nemst  glower,  radiation  chopper  and  collimating  optics,  u  is 
mounted  to  the  absorption  cell  and  evacuated  during  operation.  The  hot  gas  sample  not  only 
absorbs  but  radiates.  Locating  the  chopper  in  the  source  chamber  in  the  optical  path  before  the 
sample  cell  prevents  the  emission  from  the  hot  gas  from  being  modulated  and  recorded. 

Instrument  control,  data  acquisition  as  well  as  interferogram  data  reducUon  (including  phase 
correction,  apodization  and  Fast  Fourier  Transformation)  and  spectral  data  analysis  (including  line 
finding,  line  assignment,  line  position  measurement  and  estimation  of  molecular  constants)  are 
accomplished  with  personal  computers  in  the  laboratory  using  computer  programs  wntten 
specifically  for  the  purpose. 


DATA 

Three  interferometric  data  sets,  corresponding  to  gas  samples  at  the  same  temperature  of  800  K 
and  pressures  of  5,  15  and  40  torr,  were  used  to  obtain  the  results  reported  here.  Data  at  sample 
pressures  of  5  and  15  torr  were  collected  at  the  highest  resolution  used  in  the  study.  But  it  was 
necessary  to  apodize  this  data  to  reduce  “ringing”  in  the  wings  of  spectral  fines,  jneldmg  an 
approximate  resolution  of  0.0065  cm"'.  At  the  sample  pressure  of  40  torr,  interferometnc  data  was 
collected  using  a  mirror  displacement  of  50  cm.  Since  ringing  in  line  wings  was  not  as  big  a  problem 
at  this  pressure,  the  interferometric  data  were  transformed  unapodized,  yielding  an  approximate 
resolution  of  0.005  cm''.  The  data  collection  procedure  was  the  following:  first  a  low  resolution 
empty-cell  background  interferogram  was  collected,  then  a  high-resolution  interferogram  wth  gas 
in  the  cell  and  finally  a  second  low-resolution  empty-cell  background  interferogram.  “ 
empty-cell  background  spectra  were  in  agreement  then  they  were  co-added  and  ratioed  with  the 
high-resolution  spectrum  to  produce  a  single  transmittance  spectrum.  Multiple  data  sets  collected 
at  the  same  experimental  and  gas  sample  conditions  were  co-added  in  the  spectral  domain  as 
transmittance  spectra. 

A  representative  transmittance  spectrum  is  shown  in  Fig.  1.  This  spectrum  was  collected  at  the 
highest  resolution  employed  in  the  study.  The  lower  portion  of  Fig.  1  ‘s  a  1^'  segment  of  the 
survey  spectrum  above  it,  including  Q  branch  lines  of  one  of  the  many  C  O2  bands  which  were 
assigned  and  analyzed.  While  only  the  Q  branch  lines  are  indicated  in  the  lower  part  of  Fig.  1, 
most  of  the  remaining  lines  have  been  assigned  and  measured.  However,  because  of  the  high  density 
of  spectral  lines  some  lines,  especially  in  the  dense  region  around  15|im,  could  not  be  assigned. 

The  difficulty  of  assigning  spectral  lines,  in  dense  spectra,  to  their  proper  transitions  was  not 
unique  to  this  study  and  was  encountered  in  earlier  studies. As  a  result,  an  interactive  data 
analysis  technique  has  been  developed,'*  based  on  a  rnodmcation  of  the  method  of  Loomis  and 
Wood.'’  The  technique  has  proved  very  useful  in  identifying  and  assigning  spectral  fines  in  very 
dense  spectra.  To  implement  the  technique  however,  some  preliminary  knowledge  of  the  molecular 
constants  for  the  band  to  be  studied  must  be  available.  In  this  study  of  spectra 
that  knowledge  was  available  in  the  HITRAN*  data  base  and  particularly  the  HITEMP  data  base 
which  has  been  compiled  by  Selby'  and  Wattson  and  Rothman*’  from  results  of  theoreucal 
calculations  using  the  technique  of  Direct  Numerical  Diagonalization,  DND,  -  of  the  molecular 
Hamiltonian. 


RESULTS 

Line  positions  were  measured  for  each  of  the  more  than  6500  assigned  spectral  fines.  To  aid  the 
determination  of  spectral  line  positions,  data  points  were  interpolated  into  the  transmittance 
spectra.  This  was  achieved  by  zero  filling  the  interferograms  by  a  factor  of  16  times  their  length. 
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5  Torr  Spectrum  of  CO^  at  BOOK 

600  640  680  720  760  800  840 


Fig.  1 .  Transmittance  spectrum  of  5  torr  of  a  naturally  occurring  sample  of  carbon  dioxide  heated  to 
800  K.  The  upper  portion  of  the  figure  is  the  complete  spectrum;  the  lower  portion  is  a  detail  of  a  I  cm"‘ 
portion  showing  some  of  the  Q  branch  lines  of  the  (IIIOI  10001)  vibrational  transition  of 


For  some  rotational-vibrational  transitions,  line  positions  were  assigned  and  measured  in  two  or 
all  three  transmittance  spectra  which  were  analyzed.  In '  those  cases  the  individual  position 
measurements  were  averaged  to  yield  a  single  final  value:  that  value  was  a  weighted  average  of 
the  three  separate  measurements,  weighted  with  the  respective  position  uncertainties.  The 
corresponding  final  value  of  the  position  uncertainty  was  computed  as  the  r.m.s.  value  of  the 
separate  estimates  of  the  uncertainty.  The  estimated  uncertainty  of  individual  line  positions  was 
based  on  an  estimate  of  the  random  noise  in  the  spectrum,  the  width  and  the  strength  of  the  line. 
If  the  line  was  blended,  then  the  width  and  a  measure  of  the  asymmetry  of  the  blended  spectral 
feature  also  affected  the  estimate  of  the  uncertainty.  Besides  their  inherent  value,  the  estimated 
position  uncertainties  were  used  as  weighting  factors  in  the  least  squares  determination  of  molecular 
constants. 

Spectral  line  positions  were  calibrated  using  P,  Q  and  R  branch  lines,  for  (/)  less  than  60,  of 
the  (02201^01101)  transition.  Reference  line  positions  were  taken  from  the  HITRAN  data  base. 
They  are  based  on  the  experimental  measurements  of  Jolma  et  al.'°  For  the  spectrum  shown  in 
Fig.  1  the  average  line  position  difference  was  found  to  be  —0.00038  cm”  ^  Assigning  this  value 
to  the  band  center  frequency  of  this  transition,  667.7517  cm”',  and  assuming  a  linear  correction 
passing  through  zero  at  zero  wavenumbers,  a  linear  correction  was  made  to  all  of  the  measured 
line  positions. 

The  measured  line  positions  were  assigned  to  the  43  vibrational  transitions  indicated  in  the 
energy  level  diagram,  Fig.  2.  The  energy  levels  are  designated  by  the  notation  (uj,  /,  r)  in 
which  (v)  indicates  a  vibrational  quantum  number,  (/)  the  vibrational  angular-momentum  quantum 
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Bands  observed 


number  and  (r)  the  sequence  number  of  a  Fermi  resonance  group.^  The  expression  for  the  energy 
levels  used  in  the  least  squares  determination  of  the  molecular  constants  {Gy,  By,  Dy,Hy)  was  the 
following: 

E,(J)  =  G.  +  B[J{J  +  1)]-Z),[y(/  +  1)]^ +//,[•/(•/  +  1)1’  (*) 

where  the  subscript  (v)  indicates  the  vibrational  state  dependence  and  (J)  is  the  rotational  angular 
momentum  quantum  number.  States  for  which  (/  >  0)  exhibit  /-type  doubling  and  two  sets  of 
energy  levels,  labeled  (e)  and  (0,  result  for  these  states.  In  the  process  of  least  squares  fitting  spectral 
lines  positions,  for  bands  with  /-type  doubling,  the  following  approximate  relations  among 
molecular  constants^-^^  were  assumed  and  employed  as  constraints  for  both  upper  and  lower  states; 
if  /  =  1  then  it  is  approximately  correct  to  require  that  G.  =  Gf;  if  /  =  2  then  in  addition  --  B^, 
further,  if  /  =  3  then  in  addition  D.  =  D(,  and  finally  if  /  =  4,  then  /f,  =  //f.  The  application  of  these 
constraints  may  be  seen  in  Table  1.  For  example,  the  transition  (14402<- 13302),  centered  at 
657.6911  cm”',  has  for  the  upper  state  all  pairs  of  (e)  and  (f)  rotational--vibrational  constants 
constrained  to  be  equal  since  /  =  4,  but  for  the  lower  state  has  ^  H;,  since  /  =  3. 
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Approximately  900  of  the  6500  measured  line  positions  had  large  position  uncertainties,  usually 
because  of  line  blending.  As  a  result,  these  lines  had  very  small  weights  and  as  a  consequence  had 
very  little  effect  in  the  determination  of  molecular  constants.  For  some  bands  there  was  an 
insufficient  number  of  spectral  lines  to  estimate  H'  or  H".  This  was  determined  quantitatively  by 
comparing  the  least-squares  fit  with  and  without  H'  and  H"  included,  and  by  considering  the 
magnitude  of  the  standard  deviation  of  H'  and  H"  as  estimated  by  the  least-squares  algorithm. 


Table  1.  Rotational-vibrational  constants  (cm“')- 


Transition 

Band 

Center 

B’ 

D*xl0’ 

H’xlO^^ 

B" 

D"X10’ 

H-xlO^^ 

Range  of  Measurement 

RMS 

Error 

xio* 

20002 

lllOlc 

5942862 

38952038 

3.90141 

39034317 

3.94543 

P(0) 

R(49) 

■■ 

lllOlf 

39128975 

3.75775 

0(42) 

■■ 

11102c 

02201e 

5972382 

39075045 

1.51491 

1.2762 

39167173 

139338 

-3.4673 

P(20) 

Q(<54)  R(68) 

n 

lll02f 

02201f 

39169510 

137767 

.1922 

39167173 

139908 

3643 

P(21) 

0(75)  R(67) 

20003 

11102c 

615.8971 

39110392 

1.76174 

-11.0714 

39074060 

1.44822 

-153401 

P(45) 

R(77) 

4 

11102f 

39168504 

131260 

-14.9594 

Q(74) 

10002 

OllOle 

618.0283 

39048215 

137126 

3.6409 

39063916 

135498 

1.9061 

P(47) 

R(83) 

4 

OllOlf 

39125457 

136073 

13600 

0(102) 

21103c 

20003 

633.0965 

39101576 

139693 

-4.4813 

39110177 

1.76861 

-Z0946 

P(60) 

0(72)  R(84) 

■■ 

21103f 

39234415 

1.71558 

-5.6512 

■■ 

11102c 

10002 

647.0621 

39074843 

130989 

2.6946 

39048573 

138681 

3.9414 

P(82) 

0(98)  R(96) 

11102f 

.39169388 

138043 

3.0562 

12202c 

11102c 

6522520 

39194285 

139067 

-7.7153 

39074426 

1.48459 

.0110 

P(69) 

0(77)  R(81) 

12202f 

11102f 

39194285 

132498 

.4271 

39168975 

135475 

3831 

P(68) 

0(86)  R(76) 

01111c 

00011 

654.8692 

38759109 

134134 

38713955 

132292 

P(63) 

0(79)  R(77) 

4 

Ollllf 

38818864 

135034 

02211c 

01111c 

6552611 

38862828 

138985 

-133147 

38758857 

130372 

-4.9906 

P(46) 

0(50)  R(62) 

02211f 

Ollllf 

38862828 

130096 

-10.6679 

38818574 

130074 

-8.4453 

P(0) 

0(63)  R(61) 

13302c 

12202c 

655.6007 

39266100 

130454 

-13724 

39194618 

1.40487 

-63414 

P(58) 

0(74)  R(74) 

13302f 

12202f 

39266100 

130454 

-.4041 

39194618 

133834 

1.6257 

P(53) 

0(75)  R(71) 

14402e 

13302c 

657.6911 

39336868 

137866 

8.4642 

39266449 

135848 

8.0177 

P(39) 

0(59)  R(53) 

14402f 

13302f 

39336868 

137866 

8.4642 

39266449 

135848 

93049 

P(42) 

0(60)  R(64) 

15502 

14402 

659.2812 

39402451 

1.44331 

39332854 

1.41409 

P(38) 

0(53)  R(50) 

6 

01101c 

00001 

6673820 

39063749 

134760 

-3020 

39021820 

132899 

-3952 

P(IOO) 

0(102)  R(118) 

6 

OllOlf 

39125348 

135623 

-3236 

02201c 

OllOle 

667.7517 

39166733 

137805 

-3.2460 

39063942 

135570 

3961 

P(85) 

0(107)  R(109) 

3 

02201f 

OllOlf 

39166733 

138324 

,4566 

39125528 

136398 

.4642 

P(92) 

0(92)  R(100) 

03301c 

02201c 

668.1143 

39237963 

1.40271 

-1.4051 

39166770 

137632 

-33065 

P(88) 

Q$l)  R(98) 

3 

03301f 

02201f 

39237963 

1.40271 

-1.0203 

39166770 

138164 

3027 

P(77) 

0(89)  R(97) 

21102c 

20002 

668.2117 

39003423 

137973 

4.0871 

38955873 

134198 

6.7684 

P(60) 

0(76)  R(66) 

4 

21102f 

39116906 

23075 

04401c 

03301c 

668.4683 

39308333 

1.42338 

-1.2637 

39237956 

139955 

-1.6632 

P(67) 

0(91)  R(91) 

3 

04401f 

03301f 

39308333 

1.42338 

-13637 

39237956 

139955 

-1.2682 

P(72) 

0(90)  R(92) 

05501 

04401 

668.8133 

39378173 

1.46006 

39308650 

1.43753 

P(61) 

0(82)  R(78) 

4 

06601 

05501 

669.1490 

39445068 

1.44365 

39376251 

1.41551 

P(54) 

0(56)  R(62) 

6 

07701 

06601 

669.4789 

39504365 

33987 

39436883 

.46312 

P(36) 

0(0)  R(51) 

8 

15501 

14401 

679.0954 

39354118 

1.43111 

39289125 

139829 

P(0) 

0(50)  R(56) 

9 

14401c 

13301c 

680.0528 

39288436 

1.44956 

6.8164 

39222547 

1.42291 

4.9304 

P(49) 

0(55)  R(61) 

5 

14401f 

13301f 

39288436 

1.44956 

6.8164 

39222547 

1.42291 

73657 

P(30) 

0(62)  R(70) 

13301c 

12201c 

681.4905 

39222253 

139456 

.4304 

39155483 

1.44679 

-2.6525 

P(60) 

0(74)  R(80) 

3 

13301f 

12201f 

39222253 

139456 

2.1460 

39155483 

131285 

4.4836 

P(51) 

0(77)  R(83) 

12201e 

lllOlc 

683.8688 

39154667 

1.40730 

-8.2658 

39040845 

1.24797 

-1.4756 

P(67) 

0(91)  R(89) 

4 

12201f 

lllOlf 

39154667 

137136 

-1.4948 

39133283 

1.19999 

-2.1340 

P(64) 

0(84)  R(82) 

lllOle 

10001 

688.6717 

39041007 

136330 

1.4098 

39018913 

1.15296 

2.2521 

P(86) 

0(98)  R(IOO) 

3 

lllOlf 

39133420 

131497 

3864 

22201c 

21101c 

696.6889 

39158510 

138605 

-12.8905 

39037834 

1.16217 

1.3007 

P(33) 

0(59)  R(59) 

5 

22201f 

21101f 

39158510 

1.17285 

-1.0853 

39170651 

1.06891 

-2.5659 

P(30) 

0(52)  R(58) 

21101c 

20001 

7032363 

39037627 

1.09817 

-11.9045 

39059403 

.88971 

-12.1256 

P(62) 

0(82)  R(90) 

21101f 

39170327 

.99802 

-16.4455 

20001 

lllOle 

720.2805 

39060187 

.94981 

1.1959 

39040674 

1.23741 

-1.9230 

P(85) 

R(89) 

lllOlf 

39133157 

1.19330 

-1.9150 

0(94) 

10001 

OllOle 

720.8043 

39018418 

1.13590 

.8832 

39063425 

1.33935 

-.9092 

P(87) 

R(109) 

2 

OllOlf 

39124998 

134765 

-.8204 

0(106) 

[continued , . .] 
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Table  I — continued 


Transition 


Band 

Center 

B" 

D"xl0’ 

Range  of  Measurement 

RMS 

Error 

xlO* 

738.6731 

38956442 

135196 

5.1490 

39074769 

1.49795 

.8225 

P(71) 

R(85) 

3 

39169274 

136488 

.6246 

0(84) 

739.9474 

39038953 

1.19459 

6.7988 

39155002 

1.43548 

-2.2209 

P(50) 

0(56)  R(84) 

4 

39171775 

1.10214 

3.6859 

39155002 

130024 

P(63) 

0(75)  R(81) 

741.7244 

39040751 

1.25331 

.4825 

.39166464 

1.36868 

-4.1425 

P(86) 

0(78)  R(104) 

2 

39133160 

1.20464 

-.0919 

39166464 

137454 

-3431 

P(61) 

0(93)  R(lOl) 

7543333 

39003614 

137483 

.8039 

39194336 

139154 

-83872 

P(48) 

0(80)  R(70) 

4 

39117534 

138462 

1.6278 

.39194336 

132651 

-.2543 

P(61) 

0(81)  R(61) 

757.4788 

39154668 

1.41102 

-63531 

39237759 

139579 

-1.7773 

P(69) 

0(81)  R(83) 

3 

39154668 

1.27690 

.6288 

39237759 

139579 

-1.4634 

P(64) 

0(84)  R(80) 

7703010 

39221644 

137134 

-2.0196 

.39308332 

1.43018 

-3419 

P(48) 

0(80)  R(72) 

3 

39221644 

137134 

-.4280 

39308332 

1.43018 

-3419 

P(55) 

0(61)  R(73) 

781.7403 

39286713 

1.40032 

39376953 

1.45164 

P(46) 

0(68)  R(64) 

5 

790.9890 

39003573 

135948 

-1.4203 

39111015 

1.79570 

1.7924 

P(58) 

0(70)  R(72) 

4 

39117438 

136341 

-1.6540 

791.4476 

39041013 

135683 

.6090 

39048261 

136851 

1.7648 

P(64) 

0(82)  R(86) 

3 

39133425 

1.20894 

.1439 

828.2540 

39155669 

1.47737 

4.7331 

39075336 

135961 

13.7713 

P(51) 

0(61)  R(45) 

7 

39155669 

134451 

123306 

39169830 

1.62280 

12.8528 

P{56) 

0(54)  R(62) 

8293295 

39035264 

.94544 

38952530 

1.09295 

P(36) 

Q(54)  R(68) 

10 

39168002 

.84524 

8983472 

38863918 

138641 

1.4484 

39155000 

1.43659 

-23839 

P(64) 

0(30)  R(46) 

6 

38863918 

139445 

5.4450 

39155000 

1.29672 

53826 

P(69) 

0(0)  R(53) 

917.6463 

38709682 

1.63747 

39064023 

1.43211 

P(72) 

R(20) 

6 

927.1562 

38758726 

139061 

39040612 

1.48731 

P(69) 

0(0)  R(ll) 

3 

38818983 

132414 
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P(76) 

0(0)  R(14) 
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211011 

02211c 

022111 
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11102c 

111021 
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04401c 
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111021 

20002 

12201c 

122011 

20001 
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If  the  least-squares  fit  with  H'  and  H"  included  did  not  reduce  the  r.m.s.  line  position  error  by 
at  least  25%,  and  if  the  estimated  uncertainty  <<7,^ >  and  <<;«->  of  H'  and  H"  were  comparable  to 
or  larger  than  the  magnitude  of  the  estimates  </f'>  or  <//">,  then  these  parameters  were  omitted 
from  the  fit  Hence,  several  of  the  transitions  included  in  Table  1  have  no  estimates  of  H  and  w  . 

The  molecular  constants  given  in  Table  1  were  determined  and  are  reported  here  for  the  purpose 
of  reproducing  line  positions.  For  that  reason  the  molecular  constants  were  treated  as  fitting 
parameters  in  the  least-squares  analyses  and  not  as  physical  parameters.  As  a  consequence  the 
spectral  lines  of  each  vibrational  transition  were  least-squares  fit  independently  of  spectral  lines  o 
other  vibrational  transitions,  even  though  some  of  these  transitions  shared  common  states.  An 
for  the  same  reason,  the  parameter  standard  deviations  estimated  by  the  least  squares  algonthm 
are  not  reported  in  Table  1.  Instead,  the  number  of  significant  figures  included  in  the  table  was 
determined  by  the  precision  to  which  the  molecular  constants  may  be  used  to  predict  line  posiuons. 
The  precision  chosen  was  0.0001  cm"'  corresponding  to  a  rotational  quantum  number  (•/)  equa 
to  100.  This  criterion  was  uniformly  applied  to  all  of  the  measured  bands  and  accounts  for  the  fac 
that  the  number  of  significant  figures  in  the  table  is  the  same  for  all  of  the  5  s  and  similar  y  for 
D'%  and  5f’s.  However,  the  molecular  constants  in  Table  1  should  not  be  used  to  extopo 
calculation  of  line  positions  beyond  the  measurement  range  which  is  also  given  m  Tab  e  l. 

The  analysis  of  the  "C'^O,  fundamental  (01101^00001)  is  a  special  case  in  this  study.  To 
identify  assign  and  analyze  spectral  lines  of  weaker  bands,  the  amount  of  COj  gas  used  in  the  ce 
was  lai  enough  that  the  line  centers  of  the  low  (/)  lines  of  the  fundamental  (below  J  of  approx. 
60)  are  saturated  and  their  line  positions  are  not  well  determined.  However,  the  line  centers  of  Ae 
high  (7)  lines  (above  J  of  60)  were  not  saturated  and  their  positions  could  be  well  determined, 
molecular  constants  in  Table  1  reflect  this  fact.  For  example,  the  band  center  is  not,  and  could  no 

’'^S^em^STr^nsitions  shown  in  Fig.  1  have  not  been  analyzed  previously  and  the  line  positions 
and  molecular  constants  (for  those  bands)  measured  in  this  study  are  new.  This  includes  the 
following  transitions:  (05501  ^  04401),  (06601  -  05501),  (07701  ^06601)  as  well  as 

05502  -  14402)  and  (15501  -  14401).  Lines  of  the  transition  (08801  -  07701)  were  found  to  be 
present  in  the  figures  generated  with  the  modified  Loomis-Wood  technique  described  above,  but 
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these  lines  were  very  weak  and  the  number  of  lines  was  insufficient  to  perform  an  analysis. 
Consequently,  no  results  are  included  for  the  (08801  ^  07701)  band  in  Table  1. 

Many  of  the  line  positions  measured  in  this  study  have  been  used  to  update  the  HITRAN  data 
base;  the  latest  version  is  HITRAN92.'"  Figure  3  is  an  example  for  the  (04401  ^  03301)  band.  In 
this  figure  the  solid  circles  are  the  line  position  residuals  (observed  minus  calculated)  obtained  from 
this  study;  the  calculated  line  positions  result  from  using  the  molecular  constants  given  in  Table  1. 
The  r.m.s.  value  of  the  line  position  residuals  seen  in  Fig.  3  (3  x  lO'^^cm^’)  is  also  included  in 
Table  1.  The  solid  curves  in  Fig.  3  are  the  differences  between  the  line  positions  calculated  from 
the  molecular  constants  given  in  Table  1,  and  the  line  position  values  found  in  the  three  data  bases 
respective  labeled  “HITRAN86”,^  “Wattson”^’^'’  and  “Chedin’’;25.26  “Wattson”  refers  to  the  DND 
calculations  included  in  the  HITEMP  data  base.  As  the  curve  labeled  “HITRAN86”  shows,  the 
magnitude  of  the  improvement  is  significant  for  rotational  quantum  numbers  greater  than  +40 
and  less  than  —40.  It  should  be  noted  that  the  sharp  divergence  of  the  “HITRAN86”  curve  is 
characteristic  of  results  obtained  using  polynomial  fitting  functions,  like  equation  (1),  extrapolated 
outside  of  the  range  of  the  fit.  This  applies  to  both  the  HITRAN86  line  positions  and  the  line 
positions  calculated  from  the  molecular  constants  in  Table  1.  This  emphasizes  the  importance  of 
not  using  the  molecular  constants  in  the  table  to  extrapolate  line  positions  beyond  the  range  which 
is  also  given  there.  The  same  divergence  of  the  HITRAN86  line  positions  relative  to  the  measured 
positions  was  found  for  many  other  transitions  measured  in  this  study.  These  considerations  are 
more  appropriate  to  the  HITEMP  database,  which  includes  line  parameter  data  for  very  high 
rotational-vibrational  quantum-number  lines,  than  they  are  for  the  HITRAN  database  which 
employs  an  intensity  cutoff  criterion^  which  limits  the  vibrational  and  rotational  quantum-number 
range  of  most  bands. 

The  line  positions  measured  in  this  study  and  the  theoretical  values  which  are  found  in  the 
HITEMP  data  base  are  complementary  data  sets.  Some  of  the  line  positions  measured  in  this  study 
will  be  employed  as  input  data  in  future  DND  calculations  by  Wattson  and  Rothman;  the  results 
will  be  used  to  update  the  HITRAN  and  HITEMP  data  bases.  Conversely,  line  positions  from  the 

P  and  R  Branches 


Fig.  3.  Line  position  residuals  for  the  (04401  03301)  transition  of  (•)  are  the  line  position 

residuals  (observed  minus  calculated)  obtained  from  this  study;  the  calculated  line  positions  result  from 

using  the  molecular  constants  given  in  Table  1.  ( - are  the  differences  between  the  calculated  line 

positions  from  this  study  and  the  line  position  values  found  in  the  three  data  bases;  the  HITRAN  data 
base,  the  HITEMP  data  base  and  the  data  base  of  Chedin. 
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current  version  of  the  HITEMP  data  base  have  been  very  useful  in  this  smdy  for  the  purpose  of 
identifying  and  assigning  spectral  lines,  especially  those  lines  corresponding  to  bands  for  which 
there  was  no  data  on  the  HITRAN^  data  base  and  lines  corresponding  to  high  rotational  quantum 
numbers  (7).  This  is  illustrated  in  Fig.  3.  For  most  bands  the  DND  algorithm  leaves  a  small  shift 
compared  to  the  measurements  and  compared  to  HITRAN  for  small  values  of  the  rotational 
quantum  number.  A  general  characteristic  of  the  DND  calculations  is  that  for  high  rotational 
quantum  numbers  the  predicted  line  positions  diverge  from  the  measured  values,  but  much  more 
slowly  than  do  the  HITRAN  values.  In  addition,  for  transitions  corresponding  to  high  vibrational 
quantum  numbers,  the  calculated  line  positions  of  Wattson  and  Rothman^^®  are  in  much^better 
agreement  with  the  measured  values  than  the  calculated  positions  of  Chedin  and  Teffo. 

SUMMARY  AND  CONCLUSIONS 

The  positions  of  6500  lines  in  the  15  /zm  region  have  been  measured.  Approximately  325 

of  these  belong  to  five  previously  unmeasured  vibrational  bands.  The  remaining  measurements 
correspond  to  previously  studied  vibrational  transitions,  but  for  many  of  these  transitions  the 
measurements  of  this  study  are  new  and  extend  the  range  of  rotational  quantum  numbers  for  which 
line  position  measurements  are  now  available.  Many  of  the  position  measurements  of  this  study 
will  be  used  to  update  the  HITRAN  and  HITEMP  data  bases.  An  interactive  data  analysis 
technique  based  on  the  method  of  Loomis  and  Wood  has  proved  very  useful  in  assigning  spectral 
lines  that  otherwise  could  not  have  been  identified.  Molecular  constants  for  43  rotational— 
vibrational  transitions  of  have  been  determined  from  the  measured  line  positions. 
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Abstract — High  temperature  absorption  measurements  were  made  for  COn  gas  in  Local 
Thermodynamic  Equilibrium  (LTE)  with  a  hot  cell  and  high  resolution  interferometer.  The 
experimental  data  were  compared  to  band-model  and  line-by-Iine  model  transmittance 
calculations  using  line  parameters  from  the  HITRAN  and  HITEMP  data  bases.  The 
line-by-line  calculations  using  HITEMP  were  in  excellent  agreement  with  experimental 
measurements,  while  the  model  calculations  using  the  HITRAN  data  underpredicted  the 
absorption  by  approx.  10%. 


INTRODUCTION 

For  many  years,  the  Department  of  Defense  (DOD)  has  supported  research  on  sensor  design  for 
disturbed  atmospheric  environments.  Because  of  the  computational  resources  required  for  such  an 
effort,  these  computer  codes  have  relied  on  a  band-model  characterization  of  the  LTE  and  non-LTE 
spectral  properties  of  atmospheric  molecules,*  At  normal  temperatures  the  atmosphere  has  several 
i.r.  transmission  windows.  At  higher  temperatures,  these  windows  begin  to  close  primarily  because 
of  the  i.r.  activity  of  CO2. 

At  higher  temperatures  the  absorption  spectrum  of  COj  becomes  significantly  populated  with 
i.r.  absorption  lines  from  higher-energy,  vibrational  states.  Experimental  efforts  to  identify  and 
measure  the  absorption  strengths  of  lines  from  these  highly-excited  absorption  and  emission  bands 
are  difficult.  Thus,  in  order  to  compute  more  reasonable  high-temperature,  band-model  parameters, 
currently  it  is  necessary  to  obtain  band-strength  and  molecular-constant  data  for  highly-excited 
transitions  from  quantum  mechanical  models. 

One  of  the  major  goals  for  the  present  research  has  been  to  validate  the  theoretically  derived 
HITEMP^  data  base  for  the  atmospheric  transmission  window  near  12  /im  at  high  temperature. 
High-temperature  gas  samples  of  carbon  dioxide  were  studied  using  the  high  resolution  interferom¬ 
eter^  and  high-temperature  absorption  cell  described  in  Ref.  4.  The  spectral  region  covered  by  these 
measurements  ranged  from  600  to  1200  cm'*,  For  the  analysis  of  the  experimental  data  line- 
parameter  data  files  were  prepared  from  band-strength  and  molecular-constant  data  obtained  from 
the  HITEMP  and  HITRAN^  data  bases.  The  CO2  line  parameter  data  in  the  HITEMP  data  base 
was  calculated  by  Wattson  and  Rothman^  using  a  Direct  Numerical  Diagonalization  (DND) 
technique.  Line  parameter  data  for  CO2  was  also  obtained  from  the  1986  version  of  HITRAN.^ 
The  line  lists  were  input  to  a  line-by-line  computer  model  which  calculated  transmission  of  i.r. 
radiation  through  the  hot  cell  for  the  given  experimental  conditions.  Comparisons  with  the 
experimental  data  provided  high  temperature  validation  of  the  DND  technique.^ 

EXPERIMENTAL  DETAILS 

A  two-meter  optical  path  difference  step  and  hold  “cat’s  eye”  interferometer,  high- temperature 
gas-sample  hot  ceil  and  Nemst-glower  infrared  source  were  used  in  this  study.*^  This  interferometer 
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has  been  used  previously  to  study  two  strong  absorption  regions  in  CO2,  the  4.3  and  the  2.8  /zm 
regions.^’*  Recently,  absoq^tion  measurements  for  the  \5  region  have  been  performed  and 
analyzed  with  the  aid  of  theoretical  information  from  the  DND  model  and  many  new  band 
identifications  have  been  made.® 

The  experimental  data  described  here  were  collected  with  a  resolution  better  than  0.1  cm“^  This 
was  sufficient  for  validation  of  the  DND  data  base  for  use  with  high-temperature,  band-model, 
computer  codes.  A  single-pass,  stainless-steel,  gas-sample  celP  1.75  m  in  length  was  housed  inside 
a  commercial  electric  furnace  which  heated  the  gas  sample  to  800  K.  The  transmission  windows 
on  the  ends  of  the  cell  were  outside  of  the  furnace  enclosure,  and  maintained  at  room  temperature. 
Consequently,  there  was  a  temperature  gradient  in  the  gas  sample. 

The  experimental  data  for  a  pressure  of  760  torr  and  cell  length  of  175  cm  are  plotted  in  Fig.  1. 
The  upper  panel  is  for  a  ceil  temperature  of  300  K  and  the  lower  panel  is  for  a  temperature  of  800  K. 
Although  the  300  K  spectrum  was  not  analyzed  for  this  work,  it  serves  to  illustrate  the  dramatic 
difference  in  the  appearance  of  the  CO2  spectrum,  that  resulted  from  raising  the  temperature  of 
the  absorption  cell  from  300  to  800  K.  The  room  temperature  spectrum  consisted  of  a  single 
experimental  measurement,  whereas  the  800  K  spectrum  consisted  of  4  experimental  measurements 
coadded  together.  The  integration  time  for  the  800  K  measurements  was  also  twice  as  long  as  that 
of  the  300  K  measurements.  Both  the  300  and  the  800  K  spectra  were  ratioed  by  the  spectrum  of 
the  empty  cell,  to  remove  the  effects  of  a  nonconstant  background. 

This  experimental  apparatus  was  designed  specifically  for  line  position  measurements,  not 
intensity  measurements.  However,  since  only  a  small  portion  of  the  optical  path  experiences  a 
temperature  gradient,  the  path  was  treated  as  isothermal  for  the  relatively  low  resolution 
transmission  calculations  that  are  presented  below.  This  assumption  was  checked  by  performing 
non-isothermal  band-model calculations,  based  on  the  results  of  two-dimensional  hydrodynamics 
code  temperature  estimates.  A  more  realistic,  but  computationally  intensive,  three-dimensional 
calculation  would  predict  a  less  rapid  drop  off  in  temperature  from  the  central  portion  of  the  cell, 
and  a  sharper  drop  near  the  windows  compared  to  two-dimensional  hydrodynamics  calculations. 
Therefore,  the  two-dimensional  calculations  give  an  upper  bound  on  the  expected  temperature 
gradient  in  the  hot  cell.  These  calculations  indicate  that  for  the  experimental  conditions  of  path 
length  (175  cm),  pressure  (760  torr)  and  heater  temperature  (800  K)  the  absolute-absorption 
spectrum  was  affected  by  less  than  1%. 


RESULTS 

Experimental  absorption  data  for  a  pressure  of  760  torr  and  a  temperature  of  800  K  were 
compared  to  theoretical  transmission  calculations  using  the  HITRAN  data  and  HITEMP  data  with 
band-model  and  line-by-line  model  computer  codes.  A  series  of  comparisons  were  made  between 
the  experimental  data  and  theoretical  calculations  both  of  which  were  convolved  with  a  5cm"‘ 
FWHM  (Full  Width  at  Half  Maximum)  triangular  filter  function  to  average  over  the  rotational 
line  structure.® 

The  theoretical  expression  of  the  transmission  curve  can  be  reduced  to  an  analytic  narrow  band 
statistical  model  if  certain  simplifying  assumptions  are  made  about  the  spacing  of  the  lines  and  the 
line  strength  distribution.  For  this  work  it  is  assumed  that  all  lines  have  Lorentz  shape  and  equal 
widths,  and  that  the  distribution  of  line  strengths  is  assumed  to  follow  an  exponential  probability 
distribution  function.  The  narrow  band  statistical  model‘^  calculations  using  the  HITEMP 
band-strength  data  and  molecular  parameters  yielded  good  agreement  with  the  experimental  data 
as  shown  in  Fig.  2.  Model  calculations  using  HITRAN  data  (not  shown)  predicted  about  10%  less 
absorption  than  the  data.  This  discrepancy  is  due  to  the  many  additional  transitions  which  are 
treated  by  the  DND  model  from  which  HITEMP  is  compiled,  which  are  not  accounted  for  in  the 
HITRAN  compilation.  The  band-model  calculations  confirm  the  appropriateness  of  applying  the 
HITEMP  data-base  for  the  prediction  of  high  temperature  absorption  and  emission  characteristics 
of  CO2.  However,  band-model  calculations  rely  on  simplifying  assumptions  regarding  the 
distribution  of  line  positions  and  strengths  within  a  spectral  interval. 

Line-by-line  calculations  were  also  performed.  In  the  line-by-line  model,  line  intensities  and 
positions  were  computed  from  the  molecular  parameters  contained  in  the  HITRAN  and  HITEMP 


27 


High  temperature  absorption  measurements  and  modeling  of  CO^ 


593 
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WAVENUMBER  (enr'') 

Fig.  L  The  experimental  data  for  a  pressure  of  760  torr  and  path  length  of  175  cm  are  plotted.  In  the  upper 
panel  the  spectrum  for  7  =  300  K  are  compared  to  the  T  =  800  K  spectrum  in  the  lower  panel. 


data  bases.  For  these  line-by-line  model  calculations,  an  average  room-temperature  (296  K) 
self-broadened  width  for  CO2  lines  of  0.07  cm“’  was  adopted.’^  Assuming  a  square-root  power-law 
dependence  on  temperature,  an  average  self-broadened  half-width  of  0.04 cm"'  was  calculated  for 
a  temperature  of  800  K.  The  line-by-line  calculations  were  performed  by  generating  the  line  shape 
for  each  line  with  a  grid  resolution  of  0.015cm"'  and  carried  out  to  2.25  cm"'  (56  Lorentz  line 
widths)  on  each  side  of  line  center.  These  line-by-line  calculations  included  lines  originating  from 
bands  with  an  integrated  band  intensity  greater  than  10”®  atm/cm\  With  this  cutoff  limit  and 
limiting  the  maximum  J  to  less  than  200,  more  than  20,000  lines  were  processed  using  the 
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Fig.  3.  760torr  5.0cm‘‘  resolution  experimental  data  are  compared  to  line-by-line  (LBL)  model  calcu¬ 
lations  using  the  HITRAN  data  set. 


line-by-line  code  for  the  HITRAN  calculations,  while  more  than  50,000  CO2  lines  were  processed 
using  the  HITEMP  data.  As  shown  in  Fig.  3,  the  line-by-line  calculation  using  the  HITRAN  line 
parameter  data  exhibits  about  a  10%  predicted  decrease  in  absorption  compared  with  the 
experimental  data,  just  as  was  found  for  the  band-model  comparisons. 

As  displayed  in  Fig.  4,  there  is  good  agreement  between  the  experimental  data  and  the 
line-by-Iine  calculations  using  the  HITEMP  data.  However,  in  some  regions  of  the  spectrtim  there 
remains  a  small  discrepancy,  especially  for  the  P  and  R  branches  of  the  laser  bands  at  960  and 
1060  cm Next  the  assumption  of  constant  line  width  as  a  function  of  J  for  all  lines  in  the 
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Fig.  4.  760  torr  5.0  cm  resolution  experimental  data  are  compared  to  line-by-Iinc  model  calculations 

using  the  HITEMP  data  set. 
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Fig.  5.  The  /-dependent  half-widlii  line-by-line  (LBL)  transmission  calculation  and  data  are  compared. 
The  lower  panel  displays  the  percent  difference  between  the  two  curves  shown  in  the  upper  panel. 


line-by-line  calculation  was  removed  and  the  calculation  was  performed  using  /-dependent 
rotational  line  widths.  The  variation  of  the  width  of  the  rotational  lines  as  a  function  of  the 
rotational  quanttim  number  for  T  =  296  K  was  obtained  from  the  HITRAN  data  base.  The 
/-dependent  line  widths  were  calculated  for  T  =  800  K  and  resulted  in  the  transmission  calculation 
shown  in  Fig.  5.  The  calculation  is  in  excellent  agreement  with  the  data  and  the  difference  is  well 
within  the  experimental  uncertainty  due  to  temperature  variations  in  the  cell  and  calibration 
uncertainties.  The  square-root  power-law  dependence  on  the  ratio  of  temperatures  for  the 
calculation  of  the  Lorentz  width  was  changed  to  a  3/4  power  law  as  suggested  by  Taine,*^  but  only 
marginal  improvement  restilted  from  this  change. 


SUMMARY 

This  work  resulted  in  the  use  of  the  HITEMP  data  to  generate  high  temperature  band-model 
parameters^  for  use  with  systems  codes  under  the  DNA  program.  The  HITEMP  band  intensity  lists 
were  validated  at  the  temperature  of  800  K  for  the  spectral  range  from  650  to  1200  cm“‘. 
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ABSTRACT 

Spectra  of  a  40  torr  CO2  sample  heated  to  800K  has  been  recorded  at  high  resolution  in  the  780  to  940 
cm"^  spectral  region  using  the  Air  Force  Geophysics  Laboratory  high  resolution  interferometer.  A  total 
of  8  different,  rotation-vibration  bands  were  identified.  The  measured  line  positions  were  used  in  a  least- 
squares-fit  to  obtain  updated  rotation-vibration  constants. 

1.  INTRODUCTION 

The  high  resolution  interferometer  was  used  in  conjunction  with  a  high  temperature  absorption  cell  to 
carry  out  a  continuing  measurement  program  of  high  temperature  COj  in  several  different  wavelength 
regions^*^.  The  spectral  region  covered  by  the  present  paper  (780  to  940  cm'^)  is  just  beyond  the  high 
wavenumber  edge  of  the  V2  COj  fundamental.  'The  absorptivity  of  CO2  is  very  low  in  this  spectral  region 
at  room  temperature,  but  raising  the  temperature  to  800K  populates  many  additional  rotation-vibration 
states,  greatly  increasing  the  absorptivity. 

2.  EXPERIMENTAL  DETAILS 

The  high  resolution  interferometer  is  a  Michelson  interferometer  that  employs  cat’s  eye  retroreflectors 
and  a  step  and  integrate  mirror  motion.  For  these  measurement,  the  maximum  optical  path  length  used 
was  1  meter,  resulting  in  an  unapodized  resolution  of  0.005  cm’\  Several  spectra  were  taken  under 
identical  experimental  conditions,  checked  for  consistency,  and  then  coadded  to  increase  signal  to  noise. 
In  addition  to  the  CO2  spectra,  the  spectrum  of  the  empty  cell  was  also  measured.  The  coadded  COj 
spectrum  was  then  ratioed  with  the  empty  cell  spectrum  to  remove  the  channel  spectrum  and  take  out 
the  background. 

The  40  torr  COj  sample  was  contained  in  a  single  pass  stainless  steel  high  temperature  absorption  cell. 
The  total  absorption  path  was  1.75  meters.  The  central  1  meter  section  of  the  cell  was  maintained  at 
800K,  while  the  temperature  dropped  to  near  room  temperature  at  the  ends  of  the  cell  where  ZnSe 
windows  were  located.  Since  only  measurements  of  spectral  line  ppsitions,  and  not  line  strengths  or 
widths,  were  the  object  of  this  stucty,  no  attempt  was  made  to  characterize  the  temperature  gradients. 

The  sample  cell  is  connected  to  the  interferometer  chamber  with  a  short  bellows,  approximately  5  cm  in 
length,  to  allow  for  expansion  of  the  hot  cell.  This  5  cm  optical  path,  the  only  portion  of  the  optical  path 
which  is  not  evacuated,  is  purged  with  dry  nitrogen. 
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3.  RESULTS 

'Ihe  rotation-vibration  bands  that  were  identified 
in  the  experimental  spectra  are  indicated  by 
arrows  on  the  COj  energy  level  diagram  of 
figure  1.  A  least-squares  fit  to  the  measured  line 
positions  resulted  in  the  molecular  constants 
^ven  in  table  I.  Each  band  was  fit  independently 
with  no  attempt  to  make  a  global  COj  energy 
level  fit.  Figures  2  and  3  show  residual  plots  for 
the  14401  ■«-  05501  band.  Fig.  2  is  for  the  P  and 
R-b  ranches  and  Fig.  3  is  for  the  Q-branch.  The 
rotation-vibration  constants  used  in  the  1986 
HTTRAN  database^  and  those  obtained  by 
Richard  Wattson^  firom  a  DND  (dirert  iiumerical 
diagonalization)  calculation  are  also  indicated  on 
theL  figures.  The  present  results  were  not 
compared  to  the  1991  HTTRAN  database®  as  the 
present  molecular  constants  have  already  been 
incorporated  into  that  database. 
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Table  1.  Rotation- Vibration  Constants  (cm'‘) 


Transition  Band  Center  B  D  xlO  H  xlO 


B"  D"xl0’  HbrlO” 


14401 

05501 

781.7399 

21102e 

21102f 

20003 

790.9886 

lllOle 

lllOlf 

10002 

791.4472 

12201e 

11102e 

828.2535 

12201f 

11102f 

21101e 

21101f 

20002 

829.5290 

02211e 

12201e 

898.5467 

02211f 

12201f 

10011 

20001 

917.6459 

Olllle 

lllOle 

927.1557 

Ollllf 

lllOlf 

.39286673 

139985 

.39003464 

1.35614 

-1.6974 

.39117332 

136036 

-1.8779 

.39040970 

1.25693 

.6953 

.39133380 

130910 

.2471 

39155578 

1.47420 

4.1508 

.39155578 

1.34139 

11.9541 

.39035243 

.94662 

.39167979 

.84649 

.38863900 

138552 

.7975 

.38863900 

1.39409 

5.4035 

.38709724 

1.65222 

.38758787 

1.58655 

.38818930 

1.31813 

.39376912 

.39110912 

1.45111 

1.79262 

1.5491 

.39048219 

1.56870 

1.8676 

.39075242 

.39169731 

.38952507 

1.55608 

1.61884 

1.09417 

13.1083 

12.0884 

39154984 

39154984 

.39064079 

.39040662 

.39133204 

1.43620 

139659 

1.44636 

1.48330 

1.17040 

-2.8764 

5.3820 
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Figure  2.  Residual  plot  for  14401f  •*-  OSSOlf.P  and  Figure  3.  Residual  plot  for  14401e  05501f 

R-Branches  Q-Branch 
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INTRODUCTION 

A  high-resolution  interferometer  and  high^temperature  absorption  cell  have  been  used  in  a  continuing  measurement 
program  of  rotation-vibration  spectra  of  the  isotopes  of  carbon  dioxide  in  different  spectral  regions.  This  paper 
covers  results  for  the  isotopes  and  with  band  centers  in  the  spectral  region  600  to  800  cm"^ 

Only  a  few  experimental  details  will  be  given  here;  additional  details  may  be  found  in  the  literature.^  The 
experimental  components  include  an  infrared  source,  a  radiation  chopper,  a  high-temperature  single-pass  gas-sample 
absorption  cell,  a  Michelson  interferometer  with  "cat’s  eye"  mirrors  employing  a  step-and-integrate  mirror  motion, 
and  two  liquid-helium-cooled  detectors.  A  carbon  dioxide  gas  sample  was  contained  in  the  1.75  meter  single-pass 
stainless-steel  sample  cell.  The  central  one  meter  section  of  the  cell  was  maintained  at  a  temperature  of  800K,  with 
temperature  gradients  in  the  end  segments  of  the  cell,  ending  with  the  cell  windows  at  room  temperature.  For  this 
study  a  stepping  mirror  displacement  of  75  cm  was  used,  giving  a  resolution  (FWHM)  of 0.006  cm‘^  with  triangular 
apodization. 

Three  high-resolution  interferograms  were  collected  under  identical  experimental  conditions.  A  single  5  Torr  sample 
of  isotopically  enhanced  carbon  dioxide,  containing  approximately  88  percent  11  percent  and 

1  percent  ^“C^^02,  was  maintained  in  the  sample  cell  at  a  temperature  of  800K  throughout  the  data  collection.  After 
ajxxiization  the  timee  interferograms  were  Fourier  transformed,  checked  for  consistency  and  then  co-added  to  yield 
a  single  composite  spectrum.  This  spectrum  was  then  ratioed  with  a  low-resolution  empty-cell  spectrum  to  produce 
a  transmittance  spectrum  which  was  then  used  for  the  determination  of  spectral  line  positions. 

RESULTS 

Spectra  of  high  temperature  gas  samples  are  much  more  dense  with  spectral  lines  than  data  collected  at  room 
temperature.  Consequently,  techniques  have  been  developed  for  the  analyses  of  dense  spectral  data.^  In  the 
preliminary  analysis  discussed  here,  more  than  4000  spectral  line  positions  associated  with  24  rotation-vibration 
bands  of  the  ^^C^^02  isotope,  and  900  spectral  lines  associated  with  4  bands  of  the  isotope  have  been 

measured.  The  vibrational  transitions  associated  with  the  ^^C^^02  isotope  are  shown  in  Fig.  1.  The  band 
identifications  and  measured  band  centers  are  given  in  Table  1,  columns  1  and  2.  Five  of  the  ^^C^^02  transitions, 
those  marked  with  a  (t),  have  been  measured  previously  by  other  researchers,  and  data  for  these  bands  have  now 
been  incorporated  in  the  HITRAN  data  base.^  Also  included  in  Table  1  are  comparisons  of  the  observed  band 
centers  with  those  from  other  sources,  including  the  HITRAN^  data  base  (column  3),  those  determined  by  Wattson 
and  Rothman^  (column  4),  and  those  calculated  by  Chedin  and  Teffo^  (column  5).  The  comparisons  are  in  the  form 
of  differences  from  the  observed  band  center  estimates  reported  in  column  2  of  Table  1.  In  some  cases  the  Observed 
and  HITRAN  band  centers  agree  quite  well  (column  3),  including  transitions  between  higher  vibrational  states  such 
as  21101  20001  (0.0014  cm*^),  while  in  other  cases  the  agreement  is  much  worse,  for  example  21103  20003 

(-0.0372  cm’^).  This  results  because  some  of  the  energy  levels  shown  m  Fig.  1  have  been  previously  determined, 
by  us  or  other  researchers,  from  transitions  in  other  spectral  regions,  and  that  information  has  already  been 
incorporated  within  the  HITRAN  data  base.  The  agreement  between  the  Observed  band  centers  and  those  of  Wattson 
(column  4)  is  much  more  uniform  and  in  many  cases  is  comparable  to  or  better  than  the  agreement  between  the 
Observed  and  HITRAN  band  centers. 
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Table  1.  Observed  band  centers  compared  to  values 
determined  by  other  researchers 


j  "C“0, 

Comparison  With  0 
Researchers 

ther 

1  B«nd 

Obsafvad 

HITRAN91 

Wattson 

Chedin 

20002 

11101 

607.9669 

-.0076 

-.0021 

.0179 

20003 

11102 

610.9881 

-.0055 

-.0104 

-.0349 

10002 

01101 

617.3494t 

-.0004 

-.0195 

.0014 

21103 

20003 

619.7866 

-.0372 

-.0096 

-.0594 

11102 

10002 

630.7100 

-.0002 

-.0172 

-.0300 

12202 

11102 

63S.1375 

-.0026 

-.0168 

-.0575 

01111 

00011 

636.7515 

.0007 

.0085 

-.0065 

13302 

12202 

637.7206 

-.0429 

-.0061 

-.0974 

21102 

20002 

644.6391 

.0058 

-.0161 

-.0039 

01101 

00001 

648.477St 

-.0005 

.0130 

-.0075 

02201 

01101 

648.7861 1 

.0009 

-.0047 

-.0039 

03301 
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Figure  1.  Observed  ^C^^Oj  transitions. 

The  band  center  and  line  position  estimates  of  Wattson  are  based  on  a  technique  for  Dirwt  Numerical 

Diagonalization  (DND)  of  an  appropriate  Hamiltonian.’  In  this  technique  a  molecular  coordmate  system,  an 

appmxtTpotlntial  energy  sirface,  .md  a  set  of  basis  functions  are  f 

Hamiltonian  is  computed  and  then  iteratively  diagonalized  and  fit  to  a  set  of  exp^mentally  e 

vibration  energy  levels  using  a  non-linear  least-squares  fitting  routine  to  adjust  the  potential  energy  p^me  . 

ptm  the  derivL  potential  function,  previously  unknown  energy  levels  md  from  these 

So^,  may  be  predicted.  The  technique  has  been  applied  to  the  isotope  mmg  “ 

StsltemLedLlier  as  part  of  our  measurement  program.  The  band 

Wattson^  and  used  in  our  current  study  were  not  determined  usmg  expenmental  data  for  ^i„vin<T^n 

very  Uttle  data  existed  for  this  isotope.  Instead  the^^lts  for  this  isotope  were  obtamed  by  employmg  an 

isotopic  shift  and  using  only  daU  for  the  principal  isotope  Oj. 

Figure  2  shows  comparisons  of  rotational  line  positions  for  the  20002  -  11102e  transition  of  The  data  are 

minted  as  line  position  differences.  The  solid  black  dots  are  the  observed  positions  mmus  the  positions  deteramed 
fit  of  the  measured  positions.  The  expression  used  in  die  fit  to  represent  the  Ime  positio^ 
was  u(J)  =  E'„.(J)  -  E\.(J)  where  and  E%.(J)  are  given  m  terms  of  the  rotational  constants  .D^  , 

by  the  following  expression: 

E^(J)  =  [J(J+ 1)1  -  Dv 


80  /  FD6-3 


Line  position  differences  calculated 
using  our  rotational  constants  are 
compared,  in  Fig.  2,  to '  the  line 
positions  in  the  HITRAN^  data  base, 
those  determined  by  Wattson,^  and 
those  calculated  by  Qiedin^,  As  the  '  ^ 
figure  shows,  the  new  measurements  ^ 
are  an  improvement  over  the  ~ 

HITRAN  values,  for  this  band,  for  y 

rotational  quantum  numbers  above  S 
40.  The  same  pattern  was  found  for  ^ 
most  of  the  other  measured  bands; 
above  some  J  value  the  differences 
grow  large  very  fast.  However,  for 
some  of  the  measured  transitions, 
significant  differences  were  found  for 
all  values  of  J.  Because  the  results  of 
Wattson  (maiiced  *Wattson"  in  Fig.  2,  Line  position  differences  for  the  P  and  R  branches  of  the 

2)  were  determined  using  only  20002  11102e  band  of  ^C^^02* 

for  the  principal  isotope  and 

not  the  isotope  under  study,  I3cl6c^,  the  comparison  in  Fig.  2  is  very  impressive.  The  result  seen  in  Fig.  2  in 
w  ch  the  agreement,  between  the  line  positions  of  Wattson  and  the  measured  line  positions,  is  better  at  high  J 
v^ues  than  the  agreement  between  the  measured  line  positions  and  the  HITRAN  line  positions,  is  a  general 
characteristic  found  for  many  of  the  measured  transitions.  It  reflects  constraints  on  the  calculations  by  Wattson  for 
the  express  purpose  of  enhancing  the  high  J  results  for  which  the  data  in  the  HITRAN  data  base  may  not  be 
adequate,  as  the  curve  in  Fig.  2  marked  ’HITRAN*  fllustrates.  Figure  2  also  shows  the  general  pattern  for  the 
comparisons  with  the  positions  calculated  by  Chedin*;  the  agreement  is  better  for  low  J  values  than  for  high  J 
values;  but  as  the  figure  suggests,  the  agreement  is  generally  not  very  good. 

The  presentation  here  is  preliminary.  Analysis  of  spectral  data  taken  at  a  pressure  of  20  Torr  must  still  be  combmed 
wia  &e  remits  from  the  analysis  of  the  5  Torr  data.  Complete  results  and  estimates  of  band  centers,  line  positions 
and  from  these  rotational  constants  will  be  presented. 
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ABSTRACT 

A  high  resolution  step-and-hold  Fourier  transform  spectrometer  has  been  used  to  observe  the  spectrum  of  COj  samples  heated 
to  800K  in  a  high  tenmerature  absorption  cell.  The  interferograms  were  recorded  ^ing  a  i^imum  optical  path 
152  cm  resulting  in  an  ap^ized  spectral  resolution  of  0.007  cm"^  The  20  torr  isotopically  enriched  COj  sampb  consist^  of  88  % 
and  11%  The  sample  was  contained  in  a  single  pass  high  ten^rature  absorption  cell  with  a  total  abso^tion 

path  iSgth  of  1.75  meters.  A  total  of  1601  lines  belonging  to  15  rotation-vibration  bands  of  and  146  lin^  l^longing  to 

two  bands  of  were  identified  in  the  observed  spectrum.  A  weighted  least-squares  fit  was  then  used  to  obtam  inqjroved 

estimates  of  the  band  centers  and  rotational  constants. 

EXPERIMENTAL  SETUP 

The  experimental  conmonents  include  an  infrared  source,  high  temperature  gas  san^le  absorption  cell  and  high  resolution 
interferometer.  The  step-and-hold  Michelson  interferometer  en^loys  "cat’s-eye"  mirrors  stoppmg  displ^ment,  for 
study  of  76  cm.  This  corresponds  to  an  unapodized  resolution  of  ^proximately  0.007  cm"  .  The  mput  and  ou^ut  beams  of  the 
cat’s-eye  mirrors  are  laterally  displaced  making  both  accessible.  Since  the  two  beams  are  conq)lement^,  two  detectors  are  used 
and  the  signals  differenced  to  yield  the  interferometric  signal  and  reduce  the  effect  of  systematic  errors,  ^ch  as  souiw 
fluctuations  and  instrumental  drifts.  The  liquid  helium  cooled  detectors  are  housed  in  the  same  dewar  to  keep  their  experimental 
conditions  the  same.  In  practice  this  arrangement  is  very  stable  over  the  twelve  to  fifteen  hour  data  collection  period  necessary 
to  collect  a  single  high  resolution  interferogram.  For  operation  in  the  fifteen  micron  spectral  region  a  pot^iim  bronude 
splitter  and  copper  dop^  germanium  detectors  were  used.  The  interferometer  chamber,  which  also  contai^  me  transfer  optiw 
and  optical  filter  wheel,  is  evacuated  during  operation.  The  single  pass  stainless  steel  gas  san^le  cell  is  <^upled  to  toe 
interferometer  chamber  with  a  purge  expansion  bellows.  The  cell  is  housed  in  a  commercial  electnc  furnace  and  for  our  stuto^ 
was  heated  to  800K.  The  single  pass  length  is  1.75  meters.  However  only  the  central  one  meter  portion  of  the  cell  was  near  toe 
furnace  temperature,  with  tenqierature  g^ents  in  the  end  portions  and  the  zinc  selenide  cell  windows  at  room  temperature,  ine 
source  chamber  contains  a  Nemst  glower,  radiation  chopper  and  collimation  optics.  It  is  mounted  ngidly  to  me  absorption  wll 
and  evacuated  during  operation.  The  hot  gas  sample  not  only  absorbs  but  radiates.  Locating  the  chopper  m  toe  ^urce  ch^ber 
in  the  optical  path  before  the  sample  cell  prevents  the  gas  emissions  from  being  modulated  and  recorded.  The  mfrar^  signals 
are  detected,  dem^ulated  and  integrated  digitally.  The  digital  system  provides  fast  settle  down  after  a  movmg  mirror  step  and 
long  data  integration  times  during  the  moving  mirror  "hold". 

MEASUREMENTS 

Several  interferograms  of  high  temperature  (800K),  low  pressure  (5  and  20  torr)  isotopically  ei^c^  carbon  dioidde  sables 
were  collected  for  analyses  of  the  vibrational,  rotational  bands  in  the  9  to  17  micron  spectr^  region.  Each  high  re^lutmn 
mterferogram  of  the  5  torr  sample  included  approximately  400,000  data  points,  sampling  the  interferogram  evep^  sixth  NeHe 
laser  ftinge.  Figure  1  is  a  portion  of  a  transmittance  spectrum.  It  is  the  ratio  of  a  high  resolution  sp^trum  and  a  low  r^lution 
empty-cell  background  spectrum.  The  loss  of  instrumental  sensitivity  is  evidenced  near  600  cm"  .  For  comp^s<m,  the  lower 
portion  of  the  same  figure  is  a  spectrum,  collected  under  similar  conditions,  using  a  naturally  occun^san^le  of  ca^n 
The  figure  illustrates  the  differences  in  the  spectral  structure  of  the  and  C  O2  isotop^.  The  lyper  portion  of  Figure 

2  is  a  detail  of  the  spectrum  shown  in  Figure  1,  showing  a  smaller  spectral  region  the  5  torr  and  a  comparable  portion 

of  a  20  torr  spectrum.  The  slightly  higher  absorber  amount  in  the  20  torr  spectrum  allows  more  spectral  Imes  to  be  observed  with 
out  significant  increase  in  line  position  measurement  error. 

RESULTS 

Table  1  lists  tbe  vibrational  transitions  and  the  measured  band  center  of  the  bands  which  were  inv^tigate  in  ^s  study.  Also 
listed  in  the  table  is  the  quantum  number  range  of  the  measured  lines  in  the  P  and  R  branches  of  each  b^d  studi^.^as  weU  ^ 
the  number  of  measured  lines  for  each  band.  A  total  of  1601  spectral  lines  were  me^r^  and  a  total  of  O  C  O 

spectral  lines  were  measured.  Rotational  constants  are  not  given  in  the  table  but  are  available  from  *e  autho^To  ch^^  me 
analysis,  line  positions  calculated  from  me  estimated  rotational  constants  were  compared  wim  me  valu^  on  HITRAW92  ,  me 
most  recent  version  of  me  HITRAN  line  parameter  data  base.  As  an  example,  Figure  3  is  a  corapansion  for  the  well  known 
"laser"  band  wim  band  center  near  913  cm'^.  In  this  figure,  me  dots  represent  me  differences  between  me  me^ured  line  positions 
and  me  line  positions  calculated  from  me  estimated  rotational  constants.  The  solid  line,  marked  HITRAN92,  re^^OTte  me 
differences  between  me  line  positions  calculated  using  our  estimated  rotational  constants  and  me  line  positions  on  me  HIIKAINVZ 
data  base.  The  new  measurements  provide  some  improvement  for  this  well-known  band,  but  greater  improvement  tor  many  ot 
the  other  bands  listed  in  the  table. 


42 


NSMITTANCE 


600  640  680  720  760  800  840 

WAVENUMBERS  (cm*t) 


STorr 


< 


750  760  770  780  790  800 

WAVENUMBERS  (cm-i) 


Figure  1  Transmittance  spectra  of  800K,  5  Torr  carbon  Figure  2  Transmittance  spectra  of  5  Torr  and  20 
l^oxide.  Torr  samples  of 


Table  I  Rovibrational  transition  investigated  in  this  study. 
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Figure  3  Comparison  of  measured  line  positions  and 
line  positions  from  HITRAN92  and  positions  calculated 
from  rotational  constants. 
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Absolute  line  intensities  and  Herman- Wallis  type  intensity  parameters  of  ^^C‘^02  are  experimentally  measured  for 
the  first  time  for  the  (40‘^l)iv  ^  (00^0)  band  at  7460  cm"‘  and  the  (40°l)i  ^  {00°0)  band  at  7921  cm"*.  The  rotationless 
band  strengths  at  296  K  are  respectively 

5Sib[IV]  =  0.0429(6)  X  10'-  cm-V(moIecuIe/cm^) 


and 


=  0.00189(6)  X  10-22  cm"V(molecule/cm2). 

In  addition,  previous  measurements  on  the  (40®l)n  ^  (00°0)  and  (40°l)m  ^  (00®0)  bands  have  been  reanalyzed  to  be 
consistent  with  definitions  of  Herman- Wallis  parameters  currently  in  use.  ©  i996  Academic  Press,  inc. 


INTRODUCTION 

Recent  observations  (1-3)  of  the  nightside  near-infrared 
spectrum  of  Venus  have  discovered  emission  windows  be¬ 
tween  4040  and  4550  cm"*  and  between  5700  and  5900 
cm"*,  as  well  as  several  smaller  ones  between  7500  and 
9400  cm"*.  Carbon  dioxide  absorption  bands  are  a  major 
source  of  near-infrared  atmospheric  opacity  on  Venus;  these 
emission  windows  occur  between  CO2  bands  where  the 
opacity  is  low  enough  to  allow  radiation  from  the  hot  surface 
to  pass  through  most  of  Venus’  atmosphere.  Pollack  et  al 
(4)  analyzed  these  emission  windows  in  an  effort  to  improve 
the  determination  of  the  composition  and  cloud  structure  of 
the  lower  atmosphere  of  Venus.  To  facilitate  part  of  this 
effort.  Giver  and  Chackerian  (5)  made  laboratory  measure¬ 
ments  of  the  intensity  and  Herman -Wallis  parameters  of 
the  very  weak  (3l'0)iv  (00°0)  band  at  4416  cm"*,  which 

is  prominent  in  Venus’  emission  window  between  4040  and 
4550  cm"*.  However,  modeling  the  emission  spectrum  (4) 
did  not  produce  a  good  fit  for  the  window  centered  at  7830 
cm'*.  The  CO2  bands  that  are  prominent  in  this  region  belong 
to  the  vibrational  sequence  4i/i  +  1/3  and  associated  hot 
bands.  In  addition,  the  veiy  weak  and  perturbed  (21*2)n 
(00*^0)  perpendicular  band  contributes  some  absorption 
around  7900  cm'*. 

For  the  purpose  of  orientation  and  later  comparisons,  the 
band  centers,  rotationless  band  strengths,  and  Herman -Wal¬ 
lis  intensity  factors  for  the  series  jzi/i  +  z/3  are  taken  from 


the  1992  HITRAN  database  (6)  and  shown  in  Table  1  for  n 
=  0  through  n  =  4.  As  is  apparent  from  Table  1,  the  intensity 
sum  of  each  +  z/3  polyad  is  about  a  factor  of  50  weaker 
than  the  next  lower  sequence.  Only  two  of  the  five  bands  of 
the  (40°1)  ^  (00^0)  pentad  are  based  on  laboratory  intensity 
measurements  (7,  8).  Intensity  parameters  for  (40®l)i  ^ 
(00°0)  at  7921  cm"*,  (40®l)iv  (00*^0)  at  7460  cm"*,  and 

(40^1  )v  (00*^0)  at  7284  cm'*  have  not  been  previously 
measured.  For  the  1992  HITRAN  database,  their  rotationless 
band  strengths  and  Herman-Wallis  parameters  were  pre¬ 
dicted  from  direct  numerical  diagonalization  (DND)  calcula¬ 
tions  (9)  which  are  dependent  upon  the  5vib  values  of  mea¬ 
sured  ground  state  bands. 

For  the  purpose  of  this  paper,  we  measured  line  intensities 
to  obtain  dipole  moments,  rotationless  band  strengths,  and 
Herman-Wallis  terms  for  the  two  pentad  bands  at  7460  and 
7921  cm'*.  These  experimental  values  can  be  used  to  refine 
the  DND  calculations  for  related,  unmeasurable  bands 
needed  to  improve  the  analysis  of  Venus’  nightside  spec¬ 
trum.  Unfortunately,  the  band  at  7284  cm"  *  is  too  overlapped 
by  residual  H2O  transitions  to  be  measured  in  these  spectra. 
We  have  also  measured  the  intensities  of  some  lines  of  the 
(21*2)ii  ^  (00°0)  perpendicular  band  at  7901  cm"*  which 
exhibits  a  strong  rotational  perturbation.  As  yet,  we  have 
not  modeled  these  intensity  measurements. 

Some  intensity  measurements  were  made  of  lines  in  the 
previously  measured  bands  at  7734  and  7594  cm"*,  to  check 
the  tabulated  intensities  of  Valero  and  Boese  (7)  and  Valero 
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TABLE  1 

Positions  and  Intensities  of  the  nui  +  Bands  of 


Band 

Origin  (cm'M 

Intfpsiiv^ 

Hcnman-Wall 

is  Values^ 

^  vib 

Ai 

Aj 

(000l)i<-(0000) 

2349.143 

916076.0 

-0.00014 

0.0 

(100l)n<-(0000) 

3612.841 

9861.0 

0.00010 

1.57  »10-5 

(l00l),4-(0000) 

3714.782 

15223.0 

0.00004 

-1.14  *10-5 

(200l)ui<~(0000) 

4853.623 

78.1 

0.0008 

2.6  xlO-5 

(200l)n^(0000) 

4977.835 

347.5 

-0.00009 

-0.47  xlO-5 

(2001)i4-(0000) 

5099.660 

109.0 

0.0 

0.0 

(300l)iv<-(0000) 

6075.980 

0.514 

0.0016 

2.5  xlO-J 

(300l),„*-(0000) 

6227.917 

4.52 

0.0009 

4.2  *10-5 

(300lh<-(0000) 

6347.852 

4.54 

0.00095 

3.8  *10-5 

(300l),<-(0000) 

6503.081 

0.586 

0.00115 

3.8  *10-5 

(400l)v^(0000) 

7283.978 

0.00401 

-0.00019 

3.71  xlO-J 

(400l)iv-t-(0000) 

7460.527 

0.0469 

-0.00014 

2.82  *10-5 

'(400l)ni^(0000) 

7593.695 

0.102 

0.0 

0.0 

(40‘>1)i,4-<0000) 

7734,448 

0.0281 

0.00025 

0.0 

(400l),<-(0000) 

7920.838 

0.00175 

0.00075 

-2.43  *10-5 

^Intensities  in  units  of  lO’^^  cm-V(inolccule/cm=)  at  296K;  from  Ref.  (6). 
'’Hcrman-WalHs  parameters  defined  by  Eq.  4. 


and  lOOO-Torr  manometers  with  digital  readout,  A  total  of 
22  spectra  from  about  7200  to  8000  cm"*  were  obtained. 

In  addition  to  the  spectra  obtained  at  Ames,  a  few  spectra 
of  COo  in  the  3800  to  8400  cm"*  region  were  obtained  at 
Kitt  Peak  wi±  the  McMath  FTS,  These  spectra  utilized  a  6- 
m  base-path  White  cell  (11);  the  absorber  abundance  was 
large  enough  on  one  of  these  spectra  to  measure  the  intensi¬ 
ties  of  some  lines  of  the  7460  cm"'  band.  In  addition,  a 
second  short  cell  placed  in  series  with  the  White  cell  was 
used  for  CO  frequency  standards.  The  experimental  condi¬ 
tions  of  all  23  spectra  are  given  in  Table  2. 

Since  the  absorption  lines  in  the  7900  cm"*  region  are 
very  weak  despite  the  long  path  lengths,  we  selected  an 
aperture  that  limited  the  resolution  to  0.045  cm"*  for  most 
of  our  Ames  spectra,  which  brought  the  signal/noise  to  ac¬ 
ceptable  levels.  Interferograms  were  coadded  for  three  hours 
and  transformed  using  Hamming  apodization.  For  the 
higher-pressure  spectra  listed  in  Table  2,  this  instrumental 
resolution  was  sufficient  to  almost  fully  resolve  the  true 
spectra.  Several  of  the  lower-pressure  Ames  spectra  were 
obtained  with  0.024  cm"*  resolution.  Some  of  these  were 
divided  by  an  empty-cell  comparison  spectrum  to  yield  trans¬ 
mission  spectra.  This  eliminated  some  of  the  room  air  water 
vapor  features  from  these  spectra.  The  Kitt  Peak  spectrum 
was  obtained  with  an  hour  integration  time  using  O.'Ol  1  cm"* 


(8).  Since  our  measurements  were  within  3%  of  their  values, 
we  did  not  make  extensive  new  measurements  of  these 
bands.  However,  we  did  refit  their  measurements  in  the  same 
manner  used  in  this  paper  which  confirmed  the  1992  HI- 
TRAN  rotationless  band  strength  for  the  (40°l)iii  ^  (00^0) 
transition  and  the  revised  value  of  Valero  and  Boese  (10) 
for  the  (40°l)ii  ^  (00^0)  band.  In  this  process  we  obtained 
Herman-Wallis  parameters  consistent  with  the  definitions 
used  for  the  1992  HITRAN  database. 

EXPERIMENTAL  DETAILS 

We  obtained  spectra  using  the  Fourier  transform  spec¬ 
trometers  (FTSs)  at  NASA  Ames  and  the  Kitt  Peak  National 
Observatory.  Most  of  our  spectra  were  obtained  with  the 
Ames  25-m  base-path  White-type  (11)  multiple-traversal  ab¬ 
sorption  cell.  Operating  in  air,  radiation  from  a  quartz-halo¬ 
gen  lamp  was  sent  through  a  BOMEM  D A3. 002  Fourier- 
transform  spectrometer  equipped  with  a  quartz  beamsplitter. 
The  recombined  beam  then  entered  the  White  cell  and  subse¬ 
quently  focused  onto  an  InSb  detector  operating  at  77  K  and 
equipped  with  a  cold  filter  to  maximize  its  sensitivity  at 
>3000  cm"'.  The  cell  was  filled  with  Coleman  grade  CO2 
with  minimum  purity  of  99.99%.  The  CO2  pressures  in  the 
White  cell  were  measured  with  MKS  Baratron  10-,  I00-, 


TABLE  2 

Experimental  Conditions 


Specirum 

Number 

Path 

Length  (m) 

Pressure 

(torr) 

Temperature 

(K) 

Resolution 

(em'b 

®1 

193 

59.7 

298.0 

0.011 

2 

407 

60.3 

294.8 

0.024 

3 

407 

98.1 

294.9 

0.024 

4 

407 

138.7 

295.0 

0.024 

5 

407 

138.7 

295.0 

0.045 

6 

1207 

138.7 

294.3 

0.024 

7 

1207 

138.7 

294.7 

0.024 

8 

1207 

338.9 

295.1 

0.045 

9 

1507 

269.6 

294.7 

0.045 

10 

1507 

339.1 

294.7 

0.045 

11 

1507 

339.1 

295.1 

0.045 

12 

1607 

9.98 

294.8 

0.045 

13 

1607 

24.58 

295.2 

0.045 

14 

1607 

49.18 

295.0 

0.045 

15 

1607 

84.8 

294.8 

0.045 

16 

1607 

119.7 

295.2 

0.025 

17 

1607 

128.0 

295.1 

0.045 

18 

1607 

170.0 

294.7 

0.045 

19 

1607 

199.9 

294.7 

0.045 

20 

1607 

230.7 

295.2 

0.045 

21 

1607 

279.8 

295.7 

0.045 

22 

1707 

279.9 

294.9 

0.045 

23 

1707 

360.2 

.  295.0 

0.045 

^Spectrum 

#1  recorded  at 

the  Kitt  Peak  solar  FTS. 
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GIVER  ET  AL. 


Wavenumber  (cm“^) 

FIG.  1.  Transmittance  spectrum  of  the  (40'’l)jv  ^  (00^0)  band  at  7460 
cm"‘  obtained  at  Ames,  using  a  total  absorption  path  of  1607  m,  a  CO^ 
pressure  of  1 70.0  Torr.  and  a  temperature  of  294.7  K.  Individual  P-  and 
/?-branch  transitions  are  identified.  The  spectral  features  that  are  not  part 
of  this  band  are  primarily  incompletely  canceled  water  lines  in  the  room 
air  path  between  the  absorption  cell  and  the  Bomem  FTS. 


resolution  and  transformed  with  the  weak  “Brault”  apodiza- 
tion  (see  Spencer  et  ai  (12)).  The  (40°l)iv  (00^0)  band  at 

7460  cm”‘  is  displayed  in  Fig.  1  and  the  (40°l)i  ^  (00°0) 
band  at  7921  cm^‘  is  displayed  in  Fig.  2. 

Line  intensities  were  determined  using  nonlinear  least- 
squares  fitting  of  the  spectra  (12).  Line  profiles  were  com¬ 
puted  using  the  laboratory  conditions  for  each  spectrum,  the 
instrumental  profile  and  the  self-broadening  coefficient  for 
each  line  as  given  by  Rothman  etaL  (6).  These  self-broaden¬ 
ing  coefficients  are  a  parameterized  fit  to  several  sets  of 
pressure  broadening  measurements  on  stronger  CO2  bands, 
which  were  assumed  to  be  independent  of  vibrational  transi¬ 
tion.  These  self-broadening  coefficients  were  revised  by  up 
to  10%  from  the  1986  HITRAN  values  (13).  Since  the  pentad 
bands  we  are  measuring  are  quite  weak,  we  did  not  attempt  to 
measure  their  self-broadening  parameters.  However,  before 
accepting  the  new  coefficients  of  Rothman  et  al  (6),  we  did 
measure  the  self-broadening  coefficients  of  lines  in  the 
(30"1)|  (00°0)  band  at  6503  cm“‘  using  two  spectra  ob¬ 

tained  with  the  Kitt  Peak  FTS.  This  stronger  band  was  cho¬ 
sen  since  the  lines  have  between  30  and  80%  central  absorp¬ 
tion,  which  is  optimal  for  self-broadening  determinations; 
the  signal-noise  was  very  high  in  this  region.  These  measured 
values  agreed  with  the  1992  HITRAN  (6)  parameterization 
to  within  our  experimental  uncertainty  of  ±3%.  Thus,  these 
measurements  support  the  revision  of  the  1986  HITRAN 
self-broadening  coefficients. 

Having  adopted  the  self-broadening  coefficients  from 
Rothman  et  al.  (6),  lines  of  the  pentad  bands  were  fitted 
individually  wherever  possible.  For  each  small  spectral  inter¬ 
val,  the  position  and  intensity  of  the  line  were  adjusted  for 
the  calculated  spectrum  until  the  sum  of  the  squares  of  the 


differences  between  the  observed  and  calculated  line  profiles 
was  minimized.  In  cases  of  overlapping  lines  in  the  spectral 
interval  it  was  necessary  to  fit  two  or  more  lines  simultane¬ 
ously.  Residual  water  vapor  lines  in  the  region  of  the  7460 
cm"*  band  also  had  to  be  included  in  this  fitting  procedure; 
only  a  few  CO2  lines  had  to  be  omitted  because  of  serious 
water  vapor  contamination.  During  the  fitting  procedure  for 
each  small  spectral  interval,  the  calculated  spectrum  included 
lines  outside  the  interval,  holding  the  parameters  of  these 
lines  fixed  at  their  approximate  values. 

Line  intensities  were  determined  in  units  of  cm" '/(mole¬ 
cule/cm”)  at  the  temperature  of  each  spectrum,  using  the 
total  measured  pressure  and  path  length  without  correction 
for  isotopic  abundance.  These  line  intensities  for  each  spec¬ 
trum  were  then  standardized  to  296  K  so  the  line  measure¬ 
ments  from  different  spectra  could  be  averaged.  The  aver¬ 
aged  line  intensities,  5obs,  are  presented  in  Tables  3  and  4 
for  each  of  the  pentad  bands  measured.  These  values  are 
consistent  with  the  HITRAN  definition  for  intensities  in  a 
gas  with  a  standard  mixture  of  its  isotopomers.  Line  intensi¬ 
ties  appropriate  for  a  gas  of  the  pure  *'C'^02  isotopomer  can 
be  obtained  by  dividing  these  tabulated  values  of  5obs  by  the 
isotopic  fraction,  /  =  0.9842. 

DATA  ANALYSIS 

The  rotationless  transition  moment  squared  and  Herman- 
Wallis  intensity  parameters  for  each  band  were  obtained 
from  the  measured  line  intensities  standardized  to  296  K  via 
the  theoretical  expression  for  the  individual  line  intensities, 

=  [U^\T^V[3hcQ,r{Ty\] 

X  {o/exp(-£"/ic/itr)}Lr|/?vibl‘F(m),  [1] 

where  the  line  intensity  is  in  units  of  cm"*/(molecule/cm~). 


1.0 
I  0-9 

m 

€  0.8 


(21‘2)ji  ~(00«0)  Band 
Q  Branch  R  Branch 

2  10  20  28 

m  nTTrnrnTTn 


(10^)1, -(00^)  Band 
13c16o2P  Branch 
46  40  30 

I  I  I  1  I  i  1  I  I 


40  30  20 

-I  I  [  I  M  I  M  t  M  I 
P  Branch 


10  4 

JJLLiJ 


2  10  20  30  40 

1  m  1 1  H  n  1 1 1  It  M  1 1 1 
R  Branch 


7910  7920  7930 

Wavenumber  (cm“h 


FIG.  2.  Transmittance  spectrum  of  the  (40"  1),  (00**0)  band  at  7921 

cm”'  obtained  at  Ames,  using  a  total  absorption  path  of  1707  m,  a  CO: 
pressure  of  279.9  Torr,  and  a  temperature  of  294.9  K.  Individual  P-  and 
y?-branch  transitions  are  identified.  Q-  and  /^-branch  lines  of  the  (21'2)n-- 
(00**0)  band  at  7901  cm"'  are  also  identified  on  this  spectrum,  but  the  P- 
branch  lines  of  this  band  are  too  weak  to  observe.  In  addition.  P30  through 
P46  of  the  band  at  7981  cm”'  are  marked  on  this  spectrum. 
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TABLE  3 

Rovibrational  Line  Intensities  for  the  (40®l)i  (00®0) 

Band  of 


Line 

<T(cm'')* 

nb 

wt*= 

Sobs" 

2AS%« 

Seal*^  (Sobs'll 

cal)/Scal  (^) 

P38 

7889.024 

4 

1.81 

0.0137 

8.2 

0.0131 

5.0 

P36 

7890.793 

12 

3.22 

0.0170 

6.2 

0.0168 

1.3 

P34 

7892.553 

10 

2.60 

0.0223 

6.9 

0.021 1 

5.4 

P32 

7894.305 

11 

10.63 

0.0225 

3.4 

0.0261 

"13.6 

P30 

7896.048 

12 

11.35 

0.0344 

3.3 

0.0316 

8.7 

P28 

7897.781 

12 

9.12 

0.0352 

3.7 

0.0375 

-6.2 

P22 

7902.920 

12 

12.0 

0.0547 

1.4 

0.0549 

-0.4 

P20 

7904.611 

12 

12.0 

0.0613 

2.1 

0.0595 

3.0 

P  18 

7906.290 

12 

12.0 

0.0636 

2.7 

0.0628 

1.2 

P  14 

7909.612 

12 

12.0 

0.0655 

2.5 

0.0641 

2.1 

P  12 

7911-255 

12 

12.0 

0.0656 

1.8 

0.0615 

6.7 

P  10 

7912.885 

9 

9.0 

0.0559 

.  2.2 

0.0564 

-0.9 

P  8 

7914.502 

12 

12.0 

0.0537 

1.6 

0.0488 

9.9 

P  6 

7916.106 

12 

11.78 

0.0374 

3.2 

0.0390 

P  4 

7917.696 

10 

3.68 

0.0258 

5.8 

0,0273 

-5.5 

R  4 

7924.690 

11 

6.04 

0.0304 

4.5 

0.0348 

-12.6 

R  6 

7926.208 

12 

9.59 

0.0428 

3.6 

0,0468 

-8.5 

R  8 

7927.712 

11 

11.0 

0.0568 

2.6 

0.0570 

•0.3 

R  10 

7929.203 

U 

11.0 

0.0665 

1.8 

0.0648 

2.7 

R  12 

7930.681 

11 

11.0 

0.0713 

1.8 

0.0702 

1.5 

R  14 

7932.147 

11 

n.o 

0.0758 

1.6 

0.0731 

3.7 

R  16 

7933.599 

11 

n.o 

0.0746 

1.6 

0.0735 

1.5 

R  18 

7935.040 

10 

10.0 

0.0703 

2.2 

0.0718 

-2.1 

R  20 

7936.468 

12 

12.0 

0.0680 

2.4 

0.0682 

-0.4 

R  22 

7937.884 

12 

12.0 

0.0606 

2.3 

0.0632 

-4,2 

R24 

7939.289 

10 

lO.O 

0.0542 

2.9 

0.0573 

-5.4 

R26 

7940.682 

11 

11.0 

0.0508 

3.3 

0.0507 

0.2 

R28 

7942.065 

8 

8.0 

0.0490 

2.1 

0,0440 

11.4 

R30 

7943.439 

11 

ll.O 

0-0346 

2.3 

0.0373 

-7.4 

R32 

7944.802 

12 

11.63 

0.0310 

3.2 

0.031 1 

-0.1 

R34 

7946.156 

10 

10.0 

0.0271 

2.5 

0.0253 

7.2 

R36 

7947.501 

10 

5.22 

0.0216 

4.8 

0.0203 

6.6 

R38 

7948.838 

9 

2.46 

0.0160 

7.1 

0.0159 

0.3 

R40 

7950.167 

8 

2.67 

0.0113 

6.8 

0.0123 

-8.0 

R42 

7951.490 

3 

0.32 

0.0093 

19.5 

0.0093 

0.3 

“Line  positions  calculated  finom  Mandin  ef  al.  (15). 

*^umbcr  of  measurements. 

^Weight  used  in  least-squares  fit  for  band  intensity  parameters. 

*^Line  intensity  in  units  of  10*^  cm  •V(molecule/cm^)  at  296K. 

®Twice  the  standard  deviation  of  the  mean  of  Sobs*  2S  percentage  of  Sobs* 


the  rovibrational  partition  function  for  is  2ur(796)  — 

286.14  (14),  the  square  of  the  rotationiess  transition  moment 
has  units  of  D%  J"  is  the  lower  state  rotational  quan¬ 
tum  number,  and  /  -  0.9842  is  the  isotopic  fraction  for 
The  lower  state  rotational  energy  values  E"(cm"‘) 
were  adopted  from  the  HITRAN  tabulations,  and  the  vacuum 
line  positions  a(cm”*)  were  calculated  from  the  rotational 
constants  of  Mandin  (15).  The  Honl-London  linestrength 
factors  Lr  are  equal  to  J"  and  {J"  +  1)  respectively  for  the 
P  and  R  branches  of  parallel  ground  state  bands.  T  is  the 
Kelvin  temperature  and  k,  h,  and  c  have  their  usual  defini¬ 
tions. 

Solving  Eq.  [1]  for  where  m  =  -J"  and  J" 

-f-  1  in  the  P  and  R  branches,  respectively,  we  defined  the 
reduced  line  intensity  as  SredM  =  I  ^vib  1 

=  [3hciO^^/^7T^(7fLr}lSrQ,r{T)cxp{E''hc/kT)}.  [2] 
Following  Rothman  et  al.  (6),  the  Herman— Wallis  factor  is 
F{m)  =  [1  +  At/w  +  A2W  +  A:^m^y  [3] 


for  the  P  and  R  branches.  To  determine  the  band  intensity 
parameters,  we  found  that  cubic  fits  to  our  data  were  no 
better  than  quadratic  fits.  Thus,  we  did  weighted  quadratic 
least-squares  fits  to  the  square  root  of  the  experimentally 
determined  reduced  intensities: 

+  Aim  +  A2/n‘]  =  [5red(^)]‘^"-  ’  W 

As  suggested  by  Malathy  Devi  et  al.  (16),  the  weight  for 
each  line  in  Tables  3  and  4  was  usually  assigned  equal  to 
the  number  of  spectra  on  which  it  was  measured.  However, 
for  some  lines  the  standard  deviation  of  the  measurements 
exceeded  a  maximum  acceptable  value  6  (expressed  in  per¬ 
cent)  chosen  for  each  band.  These  lines  were  weighted  less 
than  the  lines  with  acceptable  standard  deviations;  their 
weight  was  calculated  as  the  square  of  the  ratio  of  6  divided 


TABLE  4 

Rovibrational  Line  Intensities  for  the  (40®l)iv  ^  (00®0) 
Band  of 


Line 

a{cm'*)* 

nb 

Sobs^ 

2aS%« 

Scal*^  (Sobs'S 

cal)/Scal  (^) 

P42 

7422.649 

5 

1.15 

0.248 

62 

0.236 

4.9 

P38 

7426.733 

1 

1.0 

0J98 

0.389 

2.4 

P36 

7428.734 

9 

3.21 

0.456 

3.9 

0.486 

-6.3 

P34 

7430.709 

19 

19.0 

0.607 

1.6 

0298 

KS 

P32 

743Z658 

13 

13.0 

0.695 

1.7 

0.722 

-3.7 

P30 

7434.581 

19 

19.0 

0.866 

1.3 

0.855 

1.4 

P28 

7436.478 

11 

11.0 

0.978 

1.4 

0.993 

-12 

P26 

7438.351 

14 

14.0 

1.124 

1.8 

1.131 

-0.7 

P24 

7440.199 

6 

6.0 

1258 

1.9 

1.262 

-0.3 

P22 

7442.023 

12 

12.0 

1278 

12 

1.379 

-0.0 

P20 

7443.822 

12 

12.0 

1.492 

1.4 

1.471 

1.4 

P  18 

7445.598 

11 

11.0 

1236 

2.1 

1.531 

0.3 

P  16 

7447.349 

8 

8.0 

1226 

22 

1.552 

-1.7 

P  14 

7449.078 

10 

7.36 

1240 

2.6 

1.525 

I.O 

P  12 

7450.783 

10 

10.0 

1.435 

1.6 

1.447 

-0-8 

P  10 

7452.464 

8 

5.09 

1307 

3.1 

1.314 

-0.6 

P  8 

7454.123 

14 

14.0 

1.145 

13 

1.131 

1.2 

P  6 

7455.758 

19 

19.0 

0512 

1.0 

0.898 

12 

P  4 

7457.371 

19 

19.0 

0.637 

13 

0.624 

2.0 

P  2 

7458.960 

15 

15.0 

0.322 

12 

0.321 

0.3 

R  0 

7461.302 

4 

Z72 

0.162 

42 

0.163 

-02 

R  2 

7462.834 

17 

10.82 

0.494 

2.1 

0.484 

2.1 

R  4 

7464.343 

10 

7.60 

0.797 

2.5 

0.788 

1.2 

R  6 

7465.829 

16 

16.0 

1.026 

I.l 

1.062 

■3.4 

R  8 

7467.292 

13 

13.0 

1292 

1.3 

1.295 

-0.3 

R  10 

7468.731 

8 

8.0 

1.489 

1.6 

1.479 

0.6 

R  12 

7470.147 

10 

10.0 

1.597 

1.6 

1.610 

-0.8 

R  14 

7471.539 

10 

10.0 

1.697 

1.2 

1.686 

0.7 

R  16 

7472.907 

10 

10.0 

1.703 

1.8 

1.708 

-0.3 

R  18 

7474.251 

9 

9.0 

1.696 

2.0 

1.681 

0.9 

R  20 

7475.571 

10 

10.0 

1.588 

1.9 

1.615 

-1.6 

R22 

7476.866 

7 

7.0 

1.506 

1.7 

1211 

-0.4 

R24 

7478.136 

12 

12.0 

1.356 

2.0 

1.385 

-2.1 

R26 

7479.380 

11 

11.0 

1205 

2.1 

1.242 

-3.0 

R28 

7480.599 

13 

13.0 

1.072 

12 

1.092 

'1 .8 

R  30 

7481.792 

19 

19.0 

0.963 

1.5 

0.942 

2.3 

R32 

7482.958 

14 

14.0 

0.785 

12 

0.797 

-12 

R  34 

7484.098 

14 

14.0 

0.670 

1.4 

0.661 

1.3 

R36 

7485.210 

2 

2.0 

0253 

4.4 

0.540 

2.4 

R  38 

7486.294 

16 

11.46 

0.431 

2.1 

0.432 

-0.4 

R  40 

7487.349 

17 

17.0 

0.343 

1.4 

0.340 

0.7 

R42 

7488.376 

18 

9.87 

0.262 

22 

0.264 

-0.6 

R44 

7489.373 

8 

2.85 

0.202 

4.1 

0.201 

0.7 

R46 

7490.341 

11 

6.57 

0.153 

2.7 

0.150 

1.6 

R48 

7491.277 

5 

0.73 

0.115 

8.2 

0.1 11 

4.4 

^Linc  positions  calculated  from  Mandin  et  al.  (15). 

^Number  of  measurements. 

'^Weight  used  in  least-squares  fit  for  intensity  parameters. 

'^Line  intensity  in  units  of  10'^  cm*V(molccule/cm'^)  at  296K. 

®Twice  the  standard  deviation  of  the  mean  of  Sobs,  as  percentage  of  Sobs- 
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FIG.  3.  Weighted  quadratic  least-squares  fit  of  [5red(/^0]‘^  for  the 
(40‘U)iv  {00"0)  band  to  determine  |/?vibi.  I^vibl^i*  and  i^vibl>i2- 


by  the  standard  deviation  of  the  mean  (expressed  in  per  cent) 
found  in  Tables  3' and  4.  We  chose  6  =  5.55  for  Table  3 
and  6  =  3.48  for  Table  4. 

These  fits  of  the  data  for  each  band  to  Eq.  [4]  returned 
values  and  standard  errors  for  |  /?vib  I  and  the  product  parame¬ 
ters  l^vibl^!  and  |/?vibl^2-  Figures  3  and  4  show  the  result 
of  fitting  the  square  root  of  our  reduced  experimental  intensi¬ 
ties  vs  m.  The  error  bars  on  these  figures  are  (A5%/ 
i00)[Sred(w)]‘'^,  where  the  AS%  are  half  the  values  listed  in 
column  6  of  Tables  3  and  4.  The  intensity  parameters  ob¬ 
tained  from  these  fits  are  presented  in  Table  5. 

The  rotationless  band  strengths  5?ib,  as  used  in  the  1992 
HITRAN  (6)  tabulation,  were  evaluated  at  296  K  from  these 
values  of  the  rotationless  transition  moment  of  these  ground 
state  bands, 

5^70  =  {87^M0-'V[3/lca(7)]l^ro/|/?vibl^  [5] 

using  the  value  of  the  vibrational  partition  function  QX^9(i 
K)  =  1.0846  from  Gray  and  Young  (14),  and  the  appropriate 
vibrational  band  origin,  cto  =  7920.838  and  7460.527  cm"‘ 
respectively  for  bands  I  and  IV  of  the  pentad.  These  results 
are  included  in  Table  5. 

The  uncertainties  for  the  values-  of  |  Rvib  I  were  computed 
by  two  different  methods  and  the  larger  uncertainty  was 
adopted.  In  the  first  method,  the  standard  errors  for  the  |  Rvib  I 
values  were  produced  by  the  quadratic  fitting  procedure  for 
the  averaged  values  as  shown  in  Figs.  3  and  4.  Since 
this  method  may  overlook  uncertainties  in  the  laboratory 


conditions  of  each  spectrum,  we  also  determined  |  \  val¬ 

ues  from  fits  for  12  individual  spectra  for  the  (40°l)i 
(OO^^O)  band  and  11  spectra  for  the  (40°l)fv  ^  (00°0)  band. 
These  values  were  then  averaged  for  each  band  and  the 
standard  deviation  of  the  mean  calculated.  The  line  intensity 
measurement  errors  are  the  major  source  of  uncertainty  for 
the  weaker  band:  the  goodness  of  fit  in  Fig.  4  for  the  weak 
(40‘^l)i  (00^0)  band  is  much  less  than  that  in  Fig.  3  for 

the  (40°l)iv  (00^0)  band.  Uncertainties  in  the  laboratory 

conditions,  on  the  other  hand,  have  a  nearly  equal  effect  on 
both  bands.  As  a  result,  in  Table  5  the  uncertainty  of  iRvibI 
for  the  (40°l)i  (00^0)  band  is  twice  the  standard  error  from 

the  quadratic  fit,  while  for  the  (40°l)iv  (00°0)  band  it  is 

twice  the  standard  deviation  of  the  mean  calculated  from  the 
second  method. 

The  uncertainties  in  Table  5  for  the  quantities  |i?vibl-4j 
and  I  Rvib  U2  are  the  standard  errors  from  the  quadratic  fitting 
procedure.  The  values  of  A 1  and  A2  and  their  uncertainties 
were  then  obtained  by  dividing  |/?vibUi  and  |/?vibl^2  by 
I  /?vib  I  •  Since  the  percentage  uncertainties  of  the  1 7?vib  I  values 
are  much  smaller  than  those  of  |i?vib|Ai  and  [RviblAo,  the 
uncertainties  of  Aj  and  Ao  shown  in  Table  5  have  almost 
no  dependence  on  the  |Rvibl  uncertainties.  However,  the 
determination  of  A2  is  quite  sensitive  to  the  experimental 
temperature  uncertainties.  All  the  spectra  were  obtained  with 
the  metal  White  cells  at  room  temperature;  the  temperature 
inside  is  within  1  K  of  the  temperature  outside.  To  determine 
the  dependence  of  A2  on  the  temperature,  we  repeated  the 
quadratic  fits  of  the  data  for  both  bands  assuming  a 
systematic  temperature  difference  of  1  K.  This  resulted  in 


FIG.  4.  Weighted  quadratic  least-squares  fit  of  for  the 

(40‘’l)i  ^  (00'’0)  band  to  determine  l^vibUi.  and  |/?vib|A2. 
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TABLE  5 

Intensity  Parameters  for  the  7921  and  7460  cm 
Bands  of  '^C'^02 


Parameter 

(40®l)i<-(00°0)  Band 

(40°l)iv«-(00®0)  Band 

iRvib* 

(7,95±0.I2)xl0‘^  D 

(3.902±0.029)xl0*^  D 

IRvibiAj 

(8.3±1.8)xl0'^ 

(1.78±0.23)xl0*® 

IRviblA; 

(-3.65±.98)xl0*'° 

(11.3±1.0)xl0-‘° 

“iRvibi" 

(6.32±0.20)xl0-" 

(1.523±0.023)xl0-’  D‘ 

0.001 04±0.00023 

0.00046±0.00006 

’’a. 

(.4.6±1.3)xl0'^ 

(2.9±0.5)xl0'* 

CeO 
^  vib 

(0.:89±0.006)xl0'^‘ 

(4.29±0.06)xl0-** 

roiaiionless  transition  moment  squared  in  units  of  Dcbyc^. 
*^Thc  Aj  and  A;  Hcrman-WalUs  parameters  are  dimensionless. 

^The  units  of  arc  cm‘V(molcculc/cm‘)  at  296K. 


A2  values  for  both  bands  that  differ  from  the  values  in  Table 
5  by  0.3  X  10"^  This  additional  amount  was  then  included 
in  the  uncertainty  for  the  A.  values.  The  percentage  uncer- 
tainties  for  the  values  of  |/?vibl'  and  S“i,  in  Table  5  were  set 
equal  to  twice  those  of  the  1/?,*!  values,  as  is  appropriate 
when  squaring  a  quantity. 

The  line  positions  tabulated  in  Tables  3  and  4  were  calcu¬ 
lated  from  the  rotational  constants  given  by  Mandin  (15) 
and  agree  with  the  positions  of  those  lines  tliat  are  suffi¬ 
ciently  intense  to  be  included  in  the  HITRAN  tabulation  (6). 
The  2-0  vibrational  band  of  CO  provided  calibration 
line  positions  (17)  on  the  Kitt  Peak  spectrum  of  the 
(40°l)iv  ^  (00^0)  CO2  band  (see  Table  2).  The  CO.  line 
positions  determined  from  our  fits  of  this  rpectrum  averaged 
about  0.005  cm”‘  less  than  the  HITRAN  line  positions.  This 
difference  is  well  within  the  HITRAN  uncertainty  limits  of 
0.01  to  0.1  cm*‘  for  these  line  positions;  we  conclude  from 
our  limited  set  of  calibrated  line  position  measurements  that 
the  HITRAN  positions  are  at  least  as  accurate  as  claimed. 

DISCUSSION 

The  intensities  and  Herman-Wallis  parameters  for  the 
pentad  bands  on  the  HITRAN  listing  in  Table  1  are  a  mixture 
of  values  based  on  laboratory  measurements  and  on  DND 
calculations.  The  intensities  of  the  (40^1)n  (00^0)  and 

(40”l)ni  ^  (00^0)  bands  are  based  on  the  measurements  of 
Valero  and  Boese  (7)  and  Valero  (8);  these  two  band  intensi¬ 
ties  were  included  by  Wattson  and  Rothman  (9)  in  their 
calculation  of  the  dipole-moment  function  for  DND.  There¬ 
fore,  the  HITRAN  parameters  calculated  by  DND  for  the 


other  three  pentad  bands  have  some  dependence  on  the  inten¬ 
sities  of  the  two  previously  measured  bands. 

The  intensities  of  the  (40“l),  -  (00®0)  and  (40“l),v 
(00°0)  bands  measured  on  our  spectra  are  within  10%  of  the 
HITRAN  calculated  intensities.  The  calculated  intensities 
were  revisions  by  about  20%  of  the  values  in  the  1986  HT 
TRAN  listing  of  Rothman  (13).  Both  the  1986  and  1992 
HITRAN  databases  used  DND  calculations  for  the  band  in¬ 
tensity  parameters  for  these  two  bands.  The  differences  in 
the  two  compilations  are  due  to  a  more  complete  dipole- 
moment  model  and  more  accurate  weighting  of  observed 
data  in  the  DND  calculations  for  the  1992  HITRAN  compila¬ 
tion.  This  resulted  in  decreasing  the  calculated  intensity 
for  the  (40®l)i  ^  (00°0)  band  and  increasing  it  for  the 
(40^1)iv  (00^0)  band.  Surprisingly,  our  measured  intensi¬ 

ties  for  both  bands  lie  about  halfway  between  the  calculated 
intensities  in  these  two  tabulations. 

The  Herman-Wallis  parameters  that  we  determined  in 
Table  5  can  also  be  compared  to  the  calculated  1992 
HITRAN  values  in  Table  I  for  the  (40®l)i  ^  (00°0)  and 
(40°l),v  (00^0)  bands,  since  the  DND  calculated  parame¬ 

ters  were  listed  for  unmeasured  bands.  Our  measured  A. 
coefficient  of  (2.9  ±  0.5)  x  10"^  for  the  (40°l)iv  (00^0) 

band  is  in  excellent  agreement  with  the  DND  calculated 
value  of  2.82  X  10“^  Our  measured  A.  value  of  (-4.6  ± 
1.3)  X  10”^  for  the  weak  (40°l)i  (00°0)  band  has  much 

greater  uncertainty,  but  confirms  the  DND  calculation  that 
A.  is  negative.  The  DND  calculations  for  all  these  +  vy 
polyads  show  A.  varying  smoothly  from  a  negative  value 
for  the  first  (highest  wavenumber)  member  to  the  highest 
positive  value  for  the  final  (lowest  wavenumber)  member. 


TABLE  6 

Revised  Parameters  for  the  7734  and  7594  cm 
Bands  of 


Parameter 

(40°l)ii<-(00°0)  Band 

(40°l),n4-(00®0)  Band 

iRvibi 

(3.181±0.020)xl0’^  D 

(6.02±0.07)xl0'*  D 

IRvib'A, 

(9.8±2.7)xl0'’ 

(7.3±l0.4)xl0'’ 

IRviblAj 

(-9.0±1.3)xl0-‘® 

(-5.4±4.9)xl0'*° 

“iRvibi* 

(1.012±0.013)xl0'’ 

(3.62±0.09)xl0'’  d' 

0.0003 1±0-00008 

0.00012±0.00017 

(-2.8t0.4)xl0'^ 

(-0.9±0.8)xl0'* 

'S".b 

(2.96±0.04)xl0-^^ 

(10.38±0.26)xl0'‘* 

rotationlcss  transition  moment  squared  in  units  of  Dcbyc^. 
*^Thc  Ai  and  A;  Herman-Wallis  parameters  arc  dimensionless. 

‘^Thc  units  of  arc  cm'V(nioleculc/cm^)  at  296K. 
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TABLE  7 

Some  (40*^1),  ^  (00®0)  Pentad  Band  Intensity  Parameters  Deter¬ 
mined  in  This  Work^  Divided  by  Corresponding  Parameters  Taken 
from  Table  l’’ 


i 

S^i,(ihis)/S®^b(«rrRAN) 

A 1  (this)/ A 1  (HITRAN) 

A2(this)/A2(HITRAN) 

I 

1.08 

1.39 

1.89 

II 

1.05 

1.24 

c 

III 

1.02 

c 

c 

IV 

0.91 

-329 

1.03 

^Scc  Tables  V  and  VI. 
1>1992  HITRAN  (Ref.  6). 
®Noi  determined. 


As  for  A,,  our  measured  value  of  (1.04  ±  0,23)  X  10"^  for 
the  (40^1)i  (00°0)  band  agrees  well  with  the  calculated 

value  of  0.75  X  10‘\  but  the  DND  calculation  for  the 
(40®l)iv  (00°0)  band  expects  A,  to  be  small  and  negative, 

whereas  our  measured  value  is  Ai  =  (4.6  ±  0.6)  X  10“"^. 

Valero  and  Boese  (7)  first  fit  their  line  intensity  measure¬ 
ments  with  a  linear  least-squares  and  thus  did  not  determine 
a  quadratic  Herman -Wallis  factor  for  the  (40®l)ii  (00^0) 

band.  Arie  et  al  (18)  reanalyzed  these  measurements  using 
a  quadratic  fit,  assuming  that  A  j  =  0.  This  motivated  Valero 
and  Boese  (10)  to  reanalyze  their  measurements  using  a 
quadratic  least-squares  fit  to  determine  both  linear  and  qua¬ 
dratic  Herman -Wallis  parameters;  this  fit  increased  the  in¬ 
tensity  of  the  band  by  several  percent  compared  to  their 
original  linear  fit.  However,  the  HITRAN  parameters  for 
this  band  are  based  only  on  the  original  paper  and  did  not 
incorporate  either  of  these  subsequent  quadratic  fits;  hence, 
HITRAN  assigned  A2  =  0  for  this  band.  Valero  (8)  then 
used  a  similar  quadratic  least  squares  fit  for  his  line  intensity 
measurements  of  the  (40^1)in  ^  (00°0)  band  to  obtain  the 
band  intensity  and  both  Herman -Wallis  parameters.  His 
band  intensity  was  adopted  for  the  HITRAN  tabulations; 
however,  Rothman  et  ai  (6)  assigned  both  Ai  =  0  and  A2 
—  0  for  this  band. 

Since  the  prior  fits  for  the  line  intensities  of  these  two 
bands  all  used  expressions  for  F{m)  that  differed  slightly  from 
Eq.  [3],  we  reduced  the  measurements  of  Valero  and  Boese 
(7)  and  Valero  (8)  using  the  same  procedure  as  discussed 
above  for  our  own  measurements.  These  articles  (7, 8)  tabulate 
their  line  intensities  through  7"  =  40,  which  were  determined 
from  equivalent  width  measurements  and  reduced  to  line  in¬ 
tensities  using  the  CO2  self-broadening  coefficients  of  Yama¬ 
moto  et  al  (19).  These  self-broadening  values  agree  with  the 
recent  formulation  of  Rothman  et  al  (6)  within  6%  over  J" 
=  0  through  40.  Were  it  possible  to  redetermine  their  line 
intensities  using  the  new  values  of  self-broadening  coeffi¬ 


cients,  we  estimate  that  these  intensity  values  would  change 
by  less  than  2%,  which  is  within  the  uncertainties  of  their 
individual  line  equivalent  width  measurements. 

We  made  plots  for  the  (40°l)n  ^  (OO^^O)  and  (40°l)ni 
(CX)°0)  bands  that  are  similar  to  Figs.  3  and  4  and  determined 
|y?vibL  and  |/?vibl>i2  from  quadratic  least-squares 

fits.  The  fits  for  these  bands  were  not  weighted,  since  the 
number  of  measurements  for  each  line  were  not  available 
from  Valero  and  Boese  (7)  and  Valero  (8).  Our  plot  for  the 
(40°l)ii  (00°0)  band  had  the  same  appearance  as  the  plot 

of  the  same  data  by  Arie  et  al  (18).  We  omitted  the  i^(0) 
line  for  this  fit,  since  its  measured  intensity  was  nearly  20% 
less  than  its  calculated  value,  making  it  several  times  more 
discordant  than  any  other  line.  The  parameters  we  deter¬ 
mined  and  their  corresponding  uncertainties  are  given  in 
Table  6  and  agree  with  the  values  of  |/?vibl>  and  A2 
determined  in  the  refit  of  Valero  and  Boese  (10)  when  their 
linear  and  quadratic  fitting  parameters  are  converted  to  Her- 
man-Wallis  parameters  as  defined  by  Eq.  [3].  Our  band 
intensity  value  5vib[n]  is  5%  larger  than  the  value  in  Roth¬ 
man  et  al  (6),  which  was  based  on  the  original  linear  fit  of 
Valero  and  Boese  (7). 

The  parameters  determined  from  our  fit  of  Valero ’s  (8) 
data  for  the  (40^  Om  ^  (00^0)  band  are  also  listed  in  Table 
6.  Our  band  intensity  value  Svib[in]  is  only  2%  larger  than 
the  HITRAN  value  which  was  based  on  Valero’ s  (8)  results. 
The  Herman- Wallis  parameters  we  determined  for  this  band 
have  uncertainties  of  about  the  same  magnitude;  therefore, 
HITRAN’ s  assignment  of  zero  for  these  parameters  is  not 
unreasonable.  Despite  the  close  agreement  for  the  band  in¬ 
tensity,  our  value  of  |/?vibl  disagrees  with  Valero’ s  (8)  by 
more  than  a  factor  of  3.  We  believe  he  actually  fit  |Rvibl“ 
and  his  value  of  ^sand  correctly  from  his  line  intensity  data, 
but  miscalculated  |7?vibl  from  l/?vibl  ^- 

All  these  comparisons  between  intensity  parameters  deter¬ 
mined  in  this  paper  for  the  (40^1),-  (00^0)  pentad  bands  {i 
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=  I,  n,  ni,  IV)  and  those  given  by  the  1992  HITRAN  data¬ 
base  are  summarized  as  ratios  in  Table  7. 

Having  now  detemndned  consistent  intensities  and 
Herman-Wallis  parameters  for  the  first  four  bands  of  the 
(40°1)  ^  (00^0)  pentad,  we  note  that  the  Ao  Herman-Wallis 
parameter  changes  consistently  from  a  negative  value  for 
the  first  member  of  the  pentad  to  a  positive  value  for  the 
fourth  member.  In  particular,  we  find  that  the  value  for  the 
second  and  fourth  bands  are  of  opposite  sign  and  nearly 
equal  magnitude.  This  is  consistent  with  the  DND-calculated 
A2  parameters  for  bands  of  all  these  +  1/3  polyads,  as 
are  the  measured  A2  values  for  the  dyad  bands  given  in  Table 

1.  However,  the  1992  HITRAN  values  of  A2  for  the  triad 
and  tetrad  bands  do  not  agree  with  DND  calculations.  The 
Herman-Wallis  parameters  for  the  triad  and  tetrad  bands 
are  based  on  prior  measurements;  all  the  A2  parameters  for 
the  3i/i  +  2/3  tetrad  bands  in  Table  1  have  positive  values, 
whereas  DND  calculations  of  A2  show  systematic  changes 
for  the  triad  and  tetrad  bands  similar  to  the  trends  found  for 
the  dyad  and  pentad. 

V  To  check  the  old  measured  A2  values,  we  have  made  pre¬ 
liminary  line  intensity  measurements  of  the  (30°l)i  (00°0) 

and  (30°l)iv  ^  (00°0)  bands  on  the  same  two  spectra  obtained 
at  Kitt  Peak  on  which  we  measured  self-broadening  coeffi¬ 
cients.  Reductions  of  these  new  data  show  a  negative  value 
of  A.  for  the  (30?l)i  (OO^^O)  band  and  a  positive  value  for 

the  (30®l)iv  (00°0)  band,  in  agreement  with  measurements 

done  by  Toth  et  al  (20);  this  is  consistent  with  the  DND 
calculations  for  these  bands  and  also  consistent  with  the 
trend  of  the  measured  values  for  the  dyad  and  pentad  bands. 
The  discrepancy  found  between  our  new  measurements  and 
the  1992  HITRAN  values,  which  are  based  on  the  measure¬ 
ments  of  Suarez  and  Valero  (21)  for  the  tetrad  bands,  will 
be  investigated  fully  when  additional  spectra  are  measured. 

CONCLUSION 

Our  new  measurements  of  two  bands  and  the  refitting  of 
prior  measurements  of  two  other  bands  provide  statistically 
well-determined  band  intensities  and  Herman-Wallis  pa¬ 
rameters  for  the  first  four  of  the  pentad  bands  of  ‘^C‘^02 
between  7284  and  7921  cm“^  We  recommend  that  these 
values  be  considered  for  the  next  update  of  the  HITRAN 
database.  These  new  measurements  and  analyses  have  con¬ 
firmed  the  DND  intensity  calculations  to  within  about  10% 
for  the  two  pentad  bands  not  previously  measured.  This  is 
regarded  as  good  agreement,  since  DND  is  not  expected  to 
give  exact  intensity  parameters  for  unmeasured  ground  state, 
principal  isotopomer  bands,  as  not  enough  information  is 
available  to  accurately  determine  the  dipole-moment  surface. 
However,  our  measurements  encourage  us  to  believe  that 
the  DND  intensity  parameters  for  unmeasured  bands  of  a 


polyad  are  reasonably  accurate  when  based  on  calculations 
that  include  measured  intensities  of  other  bands  of  the  same 
polyad.  Thus,  we  expect  that  the  current  1992  HITRAN 
intensity  value  for  the  still  unmeasured  (40®l)v  (00®0) 

band  at  7284  cm“^  is  accurate  to  about  10%. 
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